(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 10.1' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[    114660,       3569]
NotebookOptionsPosition[    106206,       3275]
NotebookOutlinePosition[    106579,       3291]
CellTagsIndexPosition[    106536,       3288]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{"Factor", "[", 
  RowBox[{
   SuperscriptBox["x", "2"], "-", 
   SuperscriptBox["x", 
    RowBox[{"2", " ", "n"}]], "-", 
   SuperscriptBox["x", 
    RowBox[{"4", "+", 
     RowBox[{"2", " ", "n"}]}]], "+", 
   SuperscriptBox["x", 
    RowBox[{"2", "+", 
     RowBox[{"4", " ", "n"}]}]]}], "]"}]], "Input",
 CellChangeTimes->{{3.8226983793000383`*^9, 3.8226983888106728`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(", 
   RowBox[{
    RowBox[{"-", "x"}], "+", 
    SuperscriptBox["x", "n"]}], ")"}], " ", 
  RowBox[{"(", 
   RowBox[{"x", "+", 
    SuperscriptBox["x", "n"]}], ")"}], " ", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"-", "1"}], "+", 
    SuperscriptBox["x", 
     RowBox[{"1", "+", "n"}]]}], ")"}], " ", 
  RowBox[{"(", 
   RowBox[{"1", "+", 
    SuperscriptBox["x", 
     RowBox[{"1", "+", "n"}]]}], ")"}]}]], "Output",
 CellChangeTimes->{3.822698389938961*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  SuperscriptBox["x", "2"], " ", 
  RowBox[{
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      SuperscriptBox["x", 
       RowBox[{"2", " ", "n"}]]}], ")"}], "2"], "/", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "x"}], "+", 
       SuperscriptBox["x", "n"]}], ")"}], " ", 
     RowBox[{"(", 
      RowBox[{"x", "+", 
       SuperscriptBox["x", "n"]}], ")"}], " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "1"}], "+", 
       SuperscriptBox["x", 
        RowBox[{"1", "+", "n"}]]}], ")"}], " ", 
     RowBox[{"(", 
      RowBox[{"1", "+", 
       SuperscriptBox["x", 
        RowBox[{"1", "+", "n"}]]}], ")"}]}], ")"}]}]}]], "Input",
 CellChangeTimes->{{3.822698461107071*^9, 3.8226984760434875`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SuperscriptBox["x", "2"], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      SuperscriptBox["x", 
       RowBox[{"2", " ", "n"}]]}], ")"}], "2"]}], 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "x"}], "+", 
     SuperscriptBox["x", "n"]}], ")"}], " ", 
   RowBox[{"(", 
    RowBox[{"x", "+", 
     SuperscriptBox["x", "n"]}], ")"}], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "1"}], "+", 
     SuperscriptBox["x", 
      RowBox[{"1", "+", "n"}]]}], ")"}], " ", 
   RowBox[{"(", 
    RowBox[{"1", "+", 
     SuperscriptBox["x", 
      RowBox[{"1", "+", "n"}]]}], ")"}]}]]], "Output",
 CellChangeTimes->{3.822698477144542*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FullSimplify", "[", "%", "]"}]], "Input",
 CellChangeTimes->{{3.8226984792447166`*^9, 3.8226984845083694`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SuperscriptBox["x", "2"], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      SuperscriptBox["x", 
       RowBox[{"2", " ", "n"}]]}], ")"}], "2"]}], 
  RowBox[{
   SuperscriptBox["x", "2"], "+", 
   SuperscriptBox["x", 
    RowBox[{"2", "+", 
     RowBox[{"4", " ", "n"}]}]], "-", 
   RowBox[{
    SuperscriptBox["x", 
     RowBox[{"2", " ", "n"}]], " ", 
    RowBox[{"(", 
     RowBox[{"1", "+", 
      SuperscriptBox["x", "4"]}], ")"}]}]}]]], "Output",
 CellChangeTimes->{3.82269848509476*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"eta", " ", "=", " ", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"x", "^", 
         RowBox[{"(", 
          RowBox[{"n", "+", "1"}], ")"}]}], "-", 
        RowBox[{"1", "/", 
         RowBox[{"x", "^", 
          RowBox[{"(", 
           RowBox[{"n", "+", "1"}], ")"}]}]}]}], ")"}], "/", 
      RowBox[{"(", 
       RowBox[{"x", " ", "-", " ", 
        RowBox[{"1", "/", "x"}]}], ")"}]}], ")"}], "^", "2"}], " ", "-", " ", 
   "1"}]}]], "Input",
 CellChangeTimes->{{3.822698834280723*^9, 3.822698859647937*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", "1"}], "+", 
  FractionBox[
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", 
       SuperscriptBox["x", 
        RowBox[{
         RowBox[{"-", "1"}], "-", "n"}]]}], "+", 
      SuperscriptBox["x", 
       RowBox[{"1", "+", "n"}]]}], ")"}], "2"], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", 
       FractionBox["1", "x"]}], "+", "x"}], ")"}], "2"]]}]], "Output",
 CellChangeTimes->{3.822698860943472*^9, 3.822698904136333*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData["eta"], "Input",
 CellChangeTimes->{{3.8226989052618585`*^9, 3.822698905578947*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", "1"}], "+", 
  FractionBox[
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", 
       SuperscriptBox["x", 
        RowBox[{
         RowBox[{"-", "1"}], "-", "n"}]]}], "+", 
      SuperscriptBox["x", 
       RowBox[{"1", "+", "n"}]]}], ")"}], "2"], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", 
       FractionBox["1", "x"]}], "+", "x"}], ")"}], "2"]]}]], "Output",
 CellChangeTimes->{3.8226989060631943`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"f", " ", "=", " ", 
  RowBox[{"1", "+", 
   RowBox[{"1", "/", "eta"}]}]}]], "Input",
 CellChangeTimes->{{3.8226988691977944`*^9, 3.822698875511277*^9}}],

Cell[BoxData[
 RowBox[{"1", "+", 
  FractionBox["1", 
   RowBox[{
    RowBox[{"-", "1"}], "+", 
    FractionBox[
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", 
         SuperscriptBox["x", 
          RowBox[{
           RowBox[{"-", "1"}], "-", "n"}]]}], "+", 
        SuperscriptBox["x", 
         RowBox[{"1", "+", "n"}]]}], ")"}], "2"], 
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", 
         FractionBox["1", "x"]}], "+", "x"}], ")"}], "2"]]}]]}]], "Output",
 CellChangeTimes->{3.8226988766073503`*^9, 3.822698911415497*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"f", " ", "=", " ", 
  RowBox[{"Simplify", "[", 
   RowBox[{"1", "+", 
    RowBox[{"1", "/", "eta"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8226988786075945`*^9, 3.8226988983820763`*^9}}],

Cell[BoxData[
 FractionBox[
  SuperscriptBox[
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "1"}], "+", 
     SuperscriptBox["x", 
      RowBox[{"2", "+", 
       RowBox[{"2", " ", "n"}]}]]}], ")"}], "2"], 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "1"}], "+", 
     SuperscriptBox["x", 
      RowBox[{"2", " ", "n"}]]}], ")"}], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "1"}], "+", 
     SuperscriptBox["x", 
      RowBox[{"4", "+", 
       RowBox[{"2", " ", "n"}]}]]}], ")"}]}]]], "Output",
 CellChangeTimes->{{3.822698884364526*^9, 3.8226989146081457`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"fneta", "[", "n", "]"}]], "Input",
 CellChangeTimes->{
  3.8226988907035418`*^9, {3.8226998321140003`*^9, 3.822699834769763*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", "1"}], "+", 
  RowBox[{
   SuperscriptBox[
    RowBox[{"Csch", "[", 
     FractionBox["H", 
      RowBox[{"2", " ", "T"}]], "]"}], "2"], " ", 
   SuperscriptBox[
    RowBox[{"Sinh", "[", 
     FractionBox[
      RowBox[{"H", " ", 
       RowBox[{"(", 
        RowBox[{"1", "+", "n"}], ")"}]}], 
      RowBox[{"2", " ", "T"}]], "]"}], "2"]}]}]], "Output",
 CellChangeTimes->{3.8226998350161037`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"fneta", "[", "n_", "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"Sinh", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"H", "/", "2"}], "/", "T"}], "*", 
        RowBox[{"(", 
         RowBox[{"n", "+", "1"}], ")"}]}], "]"}], "/", 
      RowBox[{"Sinh", "[", 
       RowBox[{
        RowBox[{"H", "/", "2"}], "/", "T"}], "]"}]}], ")"}], "^", "2"}], "-", 
   "1"}]}]], "Input",
 CellChangeTimes->{{3.8226996336564827`*^9, 3.822699684216403*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"an", "[", "n_", "]"}], ":=", 
  RowBox[{
   RowBox[{"Exp", "[", 
    RowBox[{
     RowBox[{"H", "/", "T"}], "*", "n"}], "]"}], "-", "1"}]}]], "Input",
 CellChangeTimes->{{3.822699691050027*^9, 3.8226996995578957`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Log", "[", 
   RowBox[{"1", "+", 
    RowBox[{"1", "/", 
     RowBox[{"fneta", "[", "n", "]"}]}]}], "]"}], "-", 
  RowBox[{"2", "*", 
   RowBox[{"Log", "[", 
    RowBox[{"an", "[", 
     RowBox[{"n", "+", "1"}], "]"}], "]"}]}], "+", 
  RowBox[{"Log", "[", 
   RowBox[{"an", "[", "n", "]"}], "]"}], "+", 
  RowBox[{"Log", "[", 
   RowBox[{"an", "[", 
    RowBox[{"n", "+", "2"}], "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8226997083101654`*^9, 3.822699755020378*^9}, {
  3.822699999424573*^9, 3.8226999995641813`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Log", "[", 
   RowBox[{
    RowBox[{"-", "1"}], "+", 
    SuperscriptBox["\[ExponentialE]", 
     FractionBox[
      RowBox[{"H", " ", "n"}], "T"]]}], "]"}], "-", 
  RowBox[{"2", " ", 
   RowBox[{"Log", "[", 
    RowBox[{
     RowBox[{"-", "1"}], "+", 
     SuperscriptBox["\[ExponentialE]", 
      FractionBox[
       RowBox[{"H", " ", 
        RowBox[{"(", 
         RowBox[{"1", "+", "n"}], ")"}]}], "T"]]}], "]"}]}], "+", 
  RowBox[{"Log", "[", 
   RowBox[{
    RowBox[{"-", "1"}], "+", 
    SuperscriptBox["\[ExponentialE]", 
     FractionBox[
      RowBox[{"H", " ", 
       RowBox[{"(", 
        RowBox[{"2", "+", "n"}], ")"}]}], "T"]]}], "]"}], "+", 
  RowBox[{"Log", "[", 
   RowBox[{"1", "+", 
    FractionBox["1", 
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      RowBox[{
       SuperscriptBox[
        RowBox[{"Csch", "[", 
         FractionBox["H", 
          RowBox[{"2", " ", "T"}]], "]"}], "2"], " ", 
       SuperscriptBox[
        RowBox[{"Sinh", "[", 
         FractionBox[
          RowBox[{"H", " ", 
           RowBox[{"(", 
            RowBox[{"1", "+", "n"}], ")"}]}], 
          RowBox[{"2", " ", "T"}]], "]"}], "2"]}]}]]}], "]"}]}]], "Output",
 CellChangeTimes->{
  3.822699756635331*^9, {3.822700000001951*^9, 3.8227000069404783`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FullSimplify", "[", "%", "]"}]], "Input",
 CellChangeTimes->{{3.822699758305174*^9, 3.8226997605835276`*^9}, {
  3.8227000093597097`*^9, 3.822700010363696*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Log", "[", 
   RowBox[{
    RowBox[{"-", "1"}], "+", 
    SuperscriptBox["\[ExponentialE]", 
     FractionBox[
      RowBox[{"H", " ", "n"}], "T"]]}], "]"}], "-", 
  RowBox[{"2", " ", 
   RowBox[{"Log", "[", 
    RowBox[{
     RowBox[{"-", "1"}], "+", 
     SuperscriptBox["\[ExponentialE]", 
      FractionBox[
       RowBox[{"H", " ", 
        RowBox[{"(", 
         RowBox[{"1", "+", "n"}], ")"}]}], "T"]]}], "]"}]}], "+", 
  RowBox[{"Log", "[", 
   RowBox[{
    RowBox[{"-", "1"}], "+", 
    SuperscriptBox["\[ExponentialE]", 
     FractionBox[
      RowBox[{"H", " ", 
       RowBox[{"(", 
        RowBox[{"2", "+", "n"}], ")"}]}], "T"]]}], "]"}], "+", 
  RowBox[{"Log", "[", 
   RowBox[{"1", "+", 
    FractionBox["1", 
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      FractionBox[
       RowBox[{
        RowBox[{"-", "1"}], "+", 
        RowBox[{"Cosh", "[", 
         FractionBox[
          RowBox[{"H", " ", 
           RowBox[{"(", 
            RowBox[{"1", "+", "n"}], ")"}]}], "T"], "]"}]}], 
       RowBox[{
        RowBox[{"-", "1"}], "+", 
        RowBox[{"Cosh", "[", 
         FractionBox["H", "T"], "]"}]}]]}]]}], "]"}]}]], "Output",
 CellChangeTimes->{
  3.8226997608647327`*^9, {3.8227000023824396`*^9, 3.8227000110927277`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"%", "/.", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"H", "\[Rule]", "1"}], ",", 
    RowBox[{"T", "\[Rule]", "2"}], ",", 
    RowBox[{"n", "\[Rule]", "4"}]}], "}"}]}]], "Input",
 CellChangeTimes->{{3.8226997679379873`*^9, 3.822699781966136*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Log", "[", 
   RowBox[{
    RowBox[{"-", "1"}], "+", 
    SuperscriptBox["\[ExponentialE]", "2"]}], "]"}], "-", 
  RowBox[{"2", " ", 
   RowBox[{"Log", "[", 
    RowBox[{
     RowBox[{"-", "1"}], "+", 
     SuperscriptBox["\[ExponentialE]", 
      RowBox[{"5", "/", "2"}]]}], "]"}]}], "+", 
  RowBox[{"Log", "[", 
   RowBox[{
    RowBox[{"-", "1"}], "+", 
    SuperscriptBox["\[ExponentialE]", "3"]}], "]"}], "+", 
  RowBox[{"Log", "[", 
   RowBox[{"1", "+", 
    FractionBox["1", 
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      FractionBox[
       RowBox[{
        RowBox[{"-", "1"}], "+", 
        RowBox[{"Cosh", "[", 
         FractionBox["5", "2"], "]"}]}], 
       RowBox[{
        RowBox[{"-", "1"}], "+", 
        RowBox[{"Cosh", "[", 
         FractionBox["1", "2"], "]"}]}]]}]]}], "]"}]}]], "Output",
 CellChangeTimes->{3.8226997824119043`*^9, 3.822700017461105*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"N", "[", "%", "]"}]], "Input",
 CellChangeTimes->{{3.82269978684406*^9, 3.8226997879941483`*^9}}],

Cell[BoxData["0.`"], "Output",
 CellChangeTimes->{3.8226997882428455`*^9, 3.8227000213299026`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"1", "+", 
  RowBox[{"1", "/", 
   RowBox[{"fneta", "[", "n", "]"}]}]}]], "Input",
 CellChangeTimes->{{3.8226998787125654`*^9, 3.822699884425299*^9}}],

Cell[BoxData[
 RowBox[{"1", "+", 
  FractionBox["1", 
   RowBox[{
    RowBox[{"-", "1"}], "+", 
    RowBox[{
     SuperscriptBox[
      RowBox[{"Csch", "[", 
       FractionBox["H", 
        RowBox[{"2", " ", "T"}]], "]"}], "2"], " ", 
     SuperscriptBox[
      RowBox[{"Sinh", "[", 
       FractionBox[
        RowBox[{"H", " ", 
         RowBox[{"(", 
          RowBox[{"1", "+", "n"}], ")"}]}], 
        RowBox[{"2", " ", "T"}]], "]"}], "2"]}]}]]}]], "Output",
 CellChangeTimes->{3.8226998864009466`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"%", "-", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"an", "[", 
      RowBox[{"n", "+", "1"}], "]"}], "^", "2"}], "/", 
    RowBox[{"an", "[", "n", "]"}]}], "/", 
   RowBox[{"an", "[", 
    RowBox[{"n", "+", "2"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.822699911222645*^9, 3.822699934053156*^9}}],

Cell[BoxData[
 RowBox[{"1", "-", 
  FractionBox[
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      SuperscriptBox["\[ExponentialE]", 
       FractionBox[
        RowBox[{"H", " ", 
         RowBox[{"(", 
          RowBox[{"1", "+", "n"}], ")"}]}], "T"]]}], ")"}], "2"], 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      SuperscriptBox["\[ExponentialE]", 
       FractionBox[
        RowBox[{"H", " ", "n"}], "T"]]}], ")"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      SuperscriptBox["\[ExponentialE]", 
       FractionBox[
        RowBox[{"H", " ", 
         RowBox[{"(", 
          RowBox[{"2", "+", "n"}], ")"}]}], "T"]]}], ")"}]}]], "+", 
  FractionBox["1", 
   RowBox[{
    RowBox[{"-", "1"}], "+", 
    RowBox[{
     SuperscriptBox[
      RowBox[{"Csch", "[", 
       FractionBox["H", 
        RowBox[{"2", " ", "T"}]], "]"}], "2"], " ", 
     SuperscriptBox[
      RowBox[{"Sinh", "[", 
       FractionBox[
        RowBox[{"H", " ", 
         RowBox[{"(", 
          RowBox[{"1", "+", "n"}], ")"}]}], 
        RowBox[{"2", " ", "T"}]], "]"}], "2"]}]}]]}]], "Output",
 CellChangeTimes->{3.8226999347930927`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Simplify", "[", "%", "]"}]], "Input",
 CellChangeTimes->{{3.8226999368740454`*^9, 3.822699939426758*^9}}],

Cell[BoxData[
 RowBox[{"1", "-", 
  FractionBox[
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      SuperscriptBox["\[ExponentialE]", 
       FractionBox[
        RowBox[{"H", " ", 
         RowBox[{"(", 
          RowBox[{"1", "+", "n"}], ")"}]}], "T"]]}], ")"}], "2"], 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      SuperscriptBox["\[ExponentialE]", 
       FractionBox[
        RowBox[{"H", " ", "n"}], "T"]]}], ")"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      SuperscriptBox["\[ExponentialE]", 
       FractionBox[
        RowBox[{"H", " ", 
         RowBox[{"(", 
          RowBox[{"2", "+", "n"}], ")"}]}], "T"]]}], ")"}]}]], "+", 
  FractionBox["1", 
   RowBox[{
    RowBox[{"-", "1"}], "+", 
    RowBox[{
     SuperscriptBox[
      RowBox[{"Csch", "[", 
       FractionBox["H", 
        RowBox[{"2", " ", "T"}]], "]"}], "2"], " ", 
     SuperscriptBox[
      RowBox[{"Sinh", "[", 
       FractionBox[
        RowBox[{"H", " ", 
         RowBox[{"(", 
          RowBox[{"1", "+", "n"}], ")"}]}], 
        RowBox[{"2", " ", "T"}]], "]"}], "2"]}]}]]}]], "Output",
 CellChangeTimes->{3.8226999401059475`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FullSimplify", "[", "%", "]"}]], "Input",
 CellChangeTimes->{{3.8226999792852783`*^9, 3.8226999834966354`*^9}}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{3.822699984038187*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"F", "[", 
   RowBox[{"x_", ",", "T_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"-", "T"}], "*", 
   RowBox[{"Log", "[", 
    RowBox[{"1", "+", 
     RowBox[{"Exp", "[", 
      RowBox[{
       RowBox[{"-", "x"}], "/", "T"}], "]"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.8227025263704815`*^9, 3.822702548412769*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"F0", "[", "x_", "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{"1", "-", 
      RowBox[{"Sign", "[", "x", "]"}]}], ")"}], "/", "2"}], "*", 
   "x"}]}]], "Input",
 CellChangeTimes->{{3.822702550084521*^9, 3.8227025739974194`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"F", "[", 
     RowBox[{"x", ",", "1"}], "]"}], "-", 
    RowBox[{"F0", "[", "x", "]"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", 
     RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8227025860651207`*^9, 3.822702603206439*^9}, {
  3.8227030032693024`*^9, 3.822703004047807*^9}}],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
     1.], LineBox[CompressedData["
1:eJxlyX1Q0wUAxvGBIC+yCKMhOmFAIITRddSa4OI0rgO5hAYVsHjVJHkX9OKl
ZAGViEQOpxBjMIFrKPJycFwO4vnB5hRwiHa8JjJ+4zeBKY0Xw7Cg9Xd/PPfc
574uiRm8z01pNJq/cf+9TwBjcWvLlqC1qriUSIiI5H3n/towWmcamycVIk/E
8Vh9ZvSujxoZN4S4rY9O1OltCUHt1/LjSiHiK2un7o4bexQz99NVIS4uew5V
tRq737pzNa8Ca/X+zb5xRlcb4uSMS5BbJqQn9b5ECEJO3uR3iVCgHtP3JdAJ
Wv/1tp3tVVhodHR4l2lDCHg2KT10CV6fCbKynbYmaHrRCXGxFA2v1K9NFloR
ArcNd4eWBij4g2UIsCRopCJn0EUGZi2TVTm7nSBYGSvehc0IeZ+f+1uJORFf
at1TqWqD84sTAq9gM0LwjmQxpagDb7TUCLKemhKs2S7rnPUuaE5LfcZkJgSL
m50V2CXHSfHFenEyjagb6T9NGnoReTbKrdl1E7T3RizTnfrhqffiR0hegJUx
MPrquhIpR8O+JbXPMW+f7nGv8A6is0O5bPafiIwdXzxUdBcVmA9D1yoiW30v
hCWMIMB04hMfk2U0HEmadNx8gOSJ/szG1ScoNuP+ugOjOFj1zS55zmO8HDyY
e9txEl9ML/8e5UAicL8mlTP0EKGa/TfuvDYFQbqqRHFUAzaqxxeLh/HQxMLv
5hUSu23MFiS+w7jEbpe21pH4JzJ1g0eqMcbbw2lqIqFcObi3O0CN8jzWSmM3
iXCPR4mlG0OoOXdcemuGRGYZa8k7cwAb1/KdJJ5aNPEbzdL4Snj/kb7V2aPF
BRmd4WKtRPzlt2R2Ki0yn53ZN/qLAsxl+nen7mnBKf/gCJehwPxw56wfqYWq
/3E5/X4fWEubzDWLOajCI9iflfWCG1iqoEXM4YAh5tijY52w7WG2Fi3NIVXe
sPClsgNXORHOjOdzkBTrM+zcO3C+Kl523YTCNsfcrwJ17SjIPGyusadQVeTp
/fH9Flio7Z6k+VMgPETlT6ca0O7qlu9UQsG921XMEtfj7wk3ASmkcD60TRYe
cxUZNdOSJjGFnw2Xi7JRC/PYOq9DbRSGK7bN7xgQoSaHe0Y0TqGyx8kkK6QC
oh+ooKxZConUgd0T6h9h/WYBm6ensJimqQxq/B7T5JW4PVsUzoIxZH+tEIdn
9labWeng6vrTRl9MPtw59QbDTh2Cb22Pi046hQ9nveI1TB0IqZo+pIhB2dsd
Cw88/u9/AbbZBmc=
      "]], LineBox[CompressedData["
1:eJwB8QMO/CFib1JlAgAAAD4AAAACAAAAhKJm57j3aD+F2OiRySHmvxuuRg9E
cnY/PBtFmtgX5r96BW3VVTSAPztso7zqDea/5rM2o4kvhT/FIwL5/wPmv1Ni
AHG9Koo/MQBfTxj65b+W30mGkhCSP4MfB0pS5uW/zLYubSyOlD+f3Uruc9zl
vwKOE1TGC5c/Cbd9rJjS5b9vPN0h+gacP54Vm3brvuW/pEy43rD+oj/HUPtA
tpflvyuD7wrAdbM/3FOiY9H85L9cOxMIrSK0P8OVboJt8uS/jPM2BZrPtD8L
pkNGDejkv+5jfv9zKbY/V56ivFfT5L+yRA30J924PzGXPmEYquS/OAYr3Y9E
vj+HYY5vSFjkv6JEs9evicQ/oJdXomK347+vx+6pf1jPPz+LmsF0gOK/vUnK
ODrC2T8A1ACQeWLgv+uqomlOWOI/uzCWeEWc3L+drS6b9LXnP8HQqpEm8ti/
vGGGFTC37D8RqNmH+OHVv3oUmdVPEvE/Y6L/d4js0r/M0NRE0przPz9N8POa
etC/ccs3Jo8W9j98sTXV+rXMv6TPRBVmyPg/06n6fzypyL+MrLeoB0z7P4hZ
f4NoW8W/AZPUScMF/j8AQDnNDD3Cvxepq8ekSABAEl3SjIttv7/Wp4AjBYgB
QBFC9LoSH7u/W6sqhnLiAkA38gs8vxW3vzybBzvFJQRAoVHjZzXXs7/jj7n2
JIQFQBZMxdkk0rC/tCN/6yHcBkCbR/J5IJWsv+CjdzIEHQhAJhLFAPaIqL/S
KEWA83gJQIcBKBs7x6S/H5pFIMi9CkCBsiPJcMihv5aqWfk5/AtACcZbgxCF
nr/Tv0LZuFUNQGL9smRa15m/a8FeCx2YDkCGbGMdIh6Wv2Kf2Za8nQ5AsDgd
JcYOlr9afVQiXKMOQFWOGrp0/5W/SDlKOZuuDkDOTIZv8eCVvySxNWcZxQ5A
aJmtQmiklb/coAzDFfIOQPscYm1FLZW/TIC6eg5MD0BkRfbWjUaUv0NeNQau
UQ9A/85MIHY4lL86PLCRTVcPQFwblRxoKpS/KPilqIxiD0DVntMTaQ6UvwRw
kdYKeQ9A/ybhRt7Wk7+9X2gyB6YPQNDoxBuQaZO/tD3jvaarD0AVhkuMEFyT
v6wbXklGsQ9A9dULSZpOk7+a11NghbwPQPAyHI7JM5O/dk8/jgPTD0CGPpKX
lv6Sv20tuhmj2A9AUx/gumDxkr9kCzWlQt4PQA59lvgz5JK/UscqvIHpD0BR
saCq9cmSv0mlpUch7w9AeIiwEuS8kr9AgyDTwPQPQJDYoHzbr5K/N2GbXmD6
D0Cellni26KSvy4/Fur//w9A267GPeWVkr/lnPI5
      "]]}, {}}, {{}, {}, {}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  Method->{
   "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None},
  PlotRange->{{-4, 4}, {-0.6922392429784324, 0.}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.8227026054179916`*^9, 3.8227030047648106`*^9, 
  3.8227272918493423`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"D", "[", 
  RowBox[{
   RowBox[{"F", "[", 
    RowBox[{"x", ",", "T"}], "]"}], ",", "T"}], "]"}]], "Input",
 CellChangeTimes->{{3.8227026249178514`*^9, 3.8227026314037375`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", 
   FractionBox[
    RowBox[{
     SuperscriptBox["\[ExponentialE]", 
      RowBox[{"-", 
       FractionBox["x", "T"]}]], " ", "x"}], 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "+", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{"-", 
         FractionBox["x", "T"]}]]}], ")"}], " ", "T"}]]}], "-", 
  RowBox[{"Log", "[", 
   RowBox[{"1", "+", 
    SuperscriptBox["\[ExponentialE]", 
     RowBox[{"-", 
      FractionBox["x", "T"]}]]}], "]"}]}]], "Output",
 CellChangeTimes->{3.822702631914816*^9, 3.8227272967492046`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Integrate", "[", 
   RowBox[{
    RowBox[{"F", "[", 
     RowBox[{"x", ",", "T"}], "]"}], ",", "x"}], "]"}], "//", 
  "FullSimplify"}]], "Input",
 CellChangeTimes->{{3.8227029449312773`*^9, 3.8227029622065687`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"-", 
    FractionBox["1", "2"]}], " ", "x", " ", 
   RowBox[{"(", 
    RowBox[{"x", "+", 
     RowBox[{"2", " ", "T", " ", 
      RowBox[{"Log", "[", 
       RowBox[{"1", "+", 
        SuperscriptBox["\[ExponentialE]", 
         RowBox[{"-", 
          FractionBox["x", "T"]}]]}], "]"}]}], "-", 
     RowBox[{"2", " ", "T", " ", 
      RowBox[{"Log", "[", 
       RowBox[{"1", "+", 
        SuperscriptBox["\[ExponentialE]", 
         FractionBox["x", "T"]]}], "]"}]}]}], ")"}]}], "+", 
  RowBox[{
   SuperscriptBox["T", "2"], " ", 
   RowBox[{"PolyLog", "[", 
    RowBox[{"2", ",", 
     RowBox[{"-", 
      SuperscriptBox["\[ExponentialE]", 
       FractionBox["x", "T"]]}]}], "]"}]}]}]], "Output",
 CellChangeTimes->{{3.8227029506641483`*^9, 3.8227029630648603`*^9}, 
   3.8227273287590065`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"dF", "[", 
   RowBox[{"x_", ",", "T_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"-", 
    FractionBox[
     RowBox[{
      SuperscriptBox["\[ExponentialE]", 
       RowBox[{"-", 
        FractionBox["x", "T"]}]], " ", "x"}], 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"1", "+", 
        SuperscriptBox["\[ExponentialE]", 
         RowBox[{"-", 
          FractionBox["x", "T"]}]]}], ")"}], " ", "T"}]]}], "-", 
   RowBox[{"Log", "[", 
    RowBox[{"1", "+", 
     SuperscriptBox["\[ExponentialE]", 
      RowBox[{"-", 
       FractionBox["x", "T"]}]]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.822702752848875*^9, 3.822702759904992*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"dF", "[", 
  RowBox[{"1", ",", "1"}], "]"}]], "Input",
 CellChangeTimes->{{3.8227027460894523`*^9, 3.8227027474728975`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", 
   FractionBox["1", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "+", 
       FractionBox["1", "\[ExponentialE]"]}], ")"}], " ", 
     "\[ExponentialE]"}]]}], "-", 
  RowBox[{"Log", "[", 
   RowBox[{"1", "+", 
    FractionBox["1", "\[ExponentialE]"]}], "]"}]}]], "Output",
 CellChangeTimes->{{3.8227027479381094`*^9, 3.822702763370681*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Integrate", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"F", "[", 
     RowBox[{"x", ",", "1"}], "]"}], "-", 
    RowBox[{"F0", "[", "x", "]"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", 
     RowBox[{"-", "140"}], ",", "140"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.822703024698206*^9, 3.822703058891302*^9}}],

Cell[BoxData[
 RowBox[{"9800", "-", 
  RowBox[{"PolyLog", "[", 
   RowBox[{"2", ",", 
    RowBox[{"-", 
     FractionBox["1", 
      SuperscriptBox["\[ExponentialE]", "140"]]}]}], "]"}], "+", 
  RowBox[{"PolyLog", "[", 
   RowBox[{"2", ",", 
    RowBox[{"-", 
     SuperscriptBox["\[ExponentialE]", "140"]}]}], "]"}]}]], "Output",
 CellChangeTimes->{{3.8227030292134404`*^9, 3.822703060661193*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"N", "[", "%", "]"}]], "Input",
 CellChangeTimes->{{3.8227030424658656`*^9, 3.8227030439415827`*^9}}],

Cell[BoxData[
 RowBox[{"-", "1.6449340668477817`"}]], "Output",
 CellChangeTimes->{{3.822703044269494*^9, 3.8227030629072657`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"NIntegrate", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"F", "[", 
      RowBox[{"x", ",", "1"}], "]"}], "-", 
     RowBox[{"F0", "[", "x", "]"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"x", ",", 
      RowBox[{"-", "Infinity"}], ",", "Infinity"}], "}"}]}], "]"}], "-", 
  RowBox[{"10000", "*", 
   RowBox[{"NIntegrate", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"F", "[", 
       RowBox[{"x", ",", "0.01"}], "]"}], "-", 
      RowBox[{"F0", "[", "x", "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", 
       RowBox[{"-", "Infinity"}], ",", "Infinity"}], "}"}]}], 
    "]"}]}]}]], "Input",
 CellChangeTimes->{{3.8227030793076534`*^9, 3.8227031085579433`*^9}, {
  3.822703195129388*^9, 3.8227032289807816`*^9}}],

Cell[BoxData[
 RowBox[{"-", "1.0652412285594437`*^-10"}]], "Output",
 CellChangeTimes->{{3.8227030992587957`*^9, 3.822703109732085*^9}, {
   3.8227031965253634`*^9, 3.8227032296937504`*^9}, 3.822727353045442*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"tf", "[", 
   RowBox[{"u_", ",", "y_"}], "]"}], ":=", 
  RowBox[{"Log", "[", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"Exp", "[", "y", "]"}], "+", 
      RowBox[{"Exp", "[", 
       RowBox[{"-", "u"}], "]"}]}], ")"}], "/", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"Exp", "[", "y", "]"}], "+", "1"}], ")"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8227054353634443`*^9, 3.822705463954508*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"tf", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "1"}], "/", 
        RowBox[{"(", 
         RowBox[{"1", "+", 
          RowBox[{"x", "^", "2"}]}], ")"}]}], ",", "0.5"}], "]"}], ",", 
     RowBox[{"1", "/", 
      RowBox[{"(", 
       RowBox[{"1", "+", 
        RowBox[{"x", "^", "2"}]}], ")"}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", 
     RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8227054668526187`*^9, 3.8227055004713745`*^9}, {
  3.822705549343641*^9, 3.8227055842495255`*^9}, {3.8227108458480263`*^9, 
  3.822710847270835*^9}, {3.8227273637342715`*^9, 3.8227273642955136`*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
    1.], LineBox[CompressedData["
1:eJw1mHc01W/8wO+1EnE1vkZJVEYpKZEG748io1AhkqxKFCUZEZEkooTMzBAy
spKs52OPjOyRez+U7GzXjN/9nfP7/fWc1znv93Pe+/2cR8T8wZXbTCQSSYlM
Iv3vKQW8YxsbFJynJ2t/j0QY6N4V91leoeDZboX01iNh4BIiLza3QMGLD+0Y
bJYLg5pxQ/OhcQreIv2Pr005DEzDY3sbuih4/kgH34ZpGATOSHyP+EzBDedc
/puKCIP5hNPpMiYUvNr7kuIZ7nAoZDe7f6eUGyfrFwQvkyLAvbFzvMyMCy84
KZ6oePY9jCYJ8J0Q3ILPHbzTeDEpGg4SapspVA7cWD64c1k0DhK3J8z3eG7G
54r3xQ/Ux0PF9frXCNhx1y5mu643CSAYKygcPsCGO2z6bvSyOgkunLvu3ObL
ikv6N2+r00mBPasWHgfUWfB0G9+KvUJpMNFxJ7HwHAu+aaTxkoJkGhRkWdZe
UGTB7R9IMpnJp8Eli7s892VY8PChspmKK2ng0WITl7ObBdcuPzK29DINfiXb
46dnmfEK8gU/sfk0+KjrSdJ6z4wbi+ajhrZ0OJwZ7WH3lwl/9PBepkpOJqiV
6FFMR5hw1rffTAfxTDBv4IrR/M2ENyYc53ndnAnhY25FEj1M+I5dSr5rE5nA
Im5Mp1Uy4aLXfLt0JT5Db6yQ9YUoJtymMiFuNu4zeAfG6u+/wIRfvtJ/aatm
FvTbx0t1ppBx8qaLI/L62fBU1kGNmkDGGx5w03GTbBCkq5kPxpDxudpN61qW
2WDgNB0y+46MO7lEb3VzzoYWZ/jH/YyM/9N9+FMkKhvKn1K/q14j40uOyulB
v7Ih0UfA6ttmMk7sj1o6YZ8DVlGBCVF3SXjbW+G1nrRcmGfX8D1xm4Sv3dJu
GM/PBQ8HpgetJiT8hvzWO2zluRCq9egUux4Jl3/cI3S1OxfKSXotj4Chb3LL
z5Q1D3be5t/Q2EHCi30t+Q6Y5UHIdcUvs082kHq34r/aPV/A4Om1fel719H8
Pc6HFkX5wInLV0byryPzvN1ZW2rzoZTMf9uXex1Zu82UFLfng+iLzo8WK//Q
pXL87qnJfJj10zko0voP+Qic5coQ+Qr+EZpHQ5/9Qwb2vF+dXn0FPE8JPPrX
ULaxfdLQzQKQGD9wXTdmFbkeGtOzli+EL16twXeDV9EtWzNCWL0Qzgo9afDw
WUXHBSrJg9cKwehyg0KG3SqK0X6jGfqkEN5+vS/MprqKhD66f0pHhbD8PPfP
16kVdK3AaylDowjqBRVsdymtIEvDG1J21sVwT+vSi1+/l9CPMWnOnF+lsK4h
3S7Zs4T8v0dNSS6WQpAqzz6HpiWUvhmYCzgRFGI/8E3fllAUryh5kyyjE2Qu
rR16s4Qkvmxc8/RB8JHv0iMn+SWkderPLXcBHH71a5tteb2Itur07tP4joPh
I20FOTk66iHpF6lrlYNfd7SZ/CEGszXtPGNYDsUKEy9O7aUjO6MD7+QtykGI
3bdJkZuOTuRYLWk/LYeBqHIT1aEFRHY++nw5sxwsq2U9DUIXUIhfdRBsqwBH
AcEaF/o8yv6kQenur4BgNHIJ5c8hURfNxxpxVRDjSuIaTZtDRnUmZRezqiDl
JH/dtvg5tOYklXsDr4LSXFUlC785JB/pUJ7WXwUjyR+PcpnOIT69gtFlkWpQ
fHtr2/XNc0iCbb+4VXI1jJr2t9GvzyJb/WzJ0ZIaAKbuq1LkGfRLSfnUZsF6
GBI3HuFZnEbmEzsnS6Tq4Y3moPPcxDTi+ys09UypHmjh01EF3dOoFt9YOnGn
Hp4e2fz7bNY0QjPeHDdy66HY6NQDfZNplMJ+tTdQ6zvY3z9JIT+ZQskvX7h3
vW2Au93ltklzE2gYF5ZjPd8MNm1mFd+oE+i7/nBLonEz2DaReJtqJpBtfFSr
vlMzOFUqFtHfT6BBWR7mpZRm8MoqZFVTnkDn04V837L9gBifnMixd+MoWlb7
3CvrH9Am/6HyiNwYCuEI8MyBFjgT8Yy/8PEwusXkmV5HbgNL6szPa3y/UNB9
HZOWhk7Q7j+UUbu/F+UE7Hlm5tALF2+Jm8JyDzKVPvehyKsXNEZEtuc39qAD
7ywPcr3rBeVpXucExx60VXn8RUhOL8iTyCpPa7tR/P7NC8enemGvcEefjHUX
8n6jedDF4icsmLhtic1tR4Sxm7fNhT6Y/e2E/vNpR1Oy998pXuuD6Tt2dv5G
7ajxVYDGhkUfjD240+3M2o7mljUytDz7gOZ+OVFXvw1VDnKfFirog9pY0TMc
qy3olehwt6owFeTQ+64xryYk31OrxDRAhZ1bWEZjZJrQCu1yjP84Ff4ZWK9c
+dWI0vQ2ash0KlTOntldBI0oMKWxMI2DBjpiNHO/le+oqcw2Z+kYDWxfC09K
2tah0i75qxfdaaDb67NBCNUhta3aQdY+NDghPsPzrrEWGRyTDbcPpME6Xiaz
drAW/cORxukEGvjPmzs3/KlGi1+VTWKqaJB6PYnF5nolSlf9JRzCSoB/Chev
CEcl4uE+eWU3FwG2Cw7iHQUVSMJKOObVDgLkA85rKPBWoJ8qliPb9xFQXT4c
wNVShrR5A3ZQFQn4RLkUX+ZehkK01gZ9VQh4bVSQ4yBVhrQ2puJ3XiRAl+7T
QfXDkUhcm86wAUNfR1fO6HUpMlJo6zjwkCH/d/+1yrEStNnomMqsI0P+5cKT
Q2olSCi4MjjalYDfhaFla0zFKJtX60rUSwJIe3suRj/+hsbEWyknogioKUp9
wNpZgP44lhzZHk/AGz2XIBuZArTzc4xKSxIBgr67uhUm89Fma9ax2UwCTk7f
uEm7mYdSeL316xAB1oWJo06VuYjTy4A9toKAGK/xB1tFc5Go+eoe/RoCmAWc
XZWHshE9PVLmfhMBcr8RE+18Nsq4R9nc2EKAZQabr1NyFjp6XHWJq4OAiOcS
knotmcjeOpsq+5OAdb78Ez5BGahXRiNHkEaAefo55SKddDR7z/7Gn34CDnUY
3xDpTEXeN4+u7hoiINBqwko3LAVFL8q0+I4QQP/n7PjSIBlp7k5V7hkjABcL
Cfjbm4gWHdtaBKcIEC3aGyUclYBOOwi088wQ8Eo7K0Xnxgd0uXleeWiWgKnf
Cl+8heJR7VHNndHzBCRPhz5/hGLRqHG2hiydkW+BnKgepWhUk6bfmLlIwKmz
jV+gMhJV2QtFbl4mgOXeSFPS+XAU7uufqrJCQFMw8whnXQhye9syb75KQHix
ENnuQjAqMK12MFtj+Pfn5M7uxrfopu5fMaV/BOgLfrJwjfNHkYF2zEzrBIzZ
9IerJb1EXhqHNiUy+Cni/b7jkyeanI+UEt8gYO/eyJWyG09Q8eNRRz8Gq1ex
mRjeeYjK9P5SfzAYj2/k+l5xA/ngYuaLDFZL2MA/SirDlXycaYPBXh8shgsW
LUGaPRANM3hX5nmKsp0jHLrBH5HNYBdWmw694x6QIun51ojBQ7lyCb5HXsC3
1MTkEYY9OuakhyUHX0FGHmlAj8E4T73ijGgAFH0ug0SGP6GEuncQTxDUXiup
aWP425AvZXPozTu4HvnIbZARD6Y323WrOcMgU4XbopMRL/nbS6dMfSOAq9v0
ZSojng/OUEVW2KJAU2Jk+MYSAR+3l7O/84qBiK31b+YY+bjv0LbI+ywOHHZ0
PL+7QEDixQJa1voHsDg86rLGyC8UHGx1U0gEq8+dAXzTBPTui67ScE0CV4kv
/DsmCeBZfZY+uJwMUilyudmjBKRZzMdmy6dCrvR035VhAs63WgQ/dfoE3I7V
nzoGCXBLveDCv5AOoR9Y7noQBPD/V2rzRyYTms/GP07pIyDPQ9osx+4zrASW
387pIWAhjWO/2FIW2KyPCVq3ESD5fIoDd8gF00m28cdVBBid3s6esTcPvor8
tjIuY8yXOTnWyOY8yGnyEBEpIWDy5tONRwfyYWp/X7JiHgE5ylzz4tQCELlT
L0+NY/T7mvTMjlffgFk4IPbzewJ2fNGdJJ8ohB81LYdNQwlwEI0a+fm2CNTZ
pFWN/BjxZZOkBiiXQjWPhXLlIwKsSrV6XWdKwSr0S4GfDQGRjnZdVjEInswV
2x+7w6i3ub4REVcc+AQGzIQMCSBu7v+92lEGPx/k+9gz5tkH5dyuz76V8NX6
BB3IBBQKq1wV66+E4nTtesoyDdrWOtuj5KogcZFILZumAcuXlRbf31UwZON8
uYGgwW3Rsw03FWrg0b65es9iGkiwteB8M3WwPc9b64MtDT7XTKa6GzQDa5VT
8NEqKtCMudyOpDbDR1lXw4eFVOCmS14ilpuh/m/G/ujPjH0SdPxJrOYPOJxn
9epbBBXe4Q9FNiZ/QMeZSQFOGypc3j1hU3q0Fe7bdscr8FChvrOfTeFrO8SN
X508r9EHJer18qfKe+CDGraLN7wXxj99sNIy+AWW1rcnVKM7ITpD6v38l2Fw
d0kb4O9h3LOF33e5bhi+uWmf7/neCor3yE7r1GGof1xTFlfaCjSJ9ivsbCPg
Sh0QVklsBeFEZw7BqyPgeJN2f/BBKyRGVjqfo4/ALas0cS/2Vkh/aWgQJDcG
XS9S2hwY+7/YzPs/6a8TcNJWdWX3q2bo4aG2/WmYgD1MT/giHjYDHckEvf81
AfmzZafErzWDtNAA9ybuvxB30y/IVaIZEnpPs9Nu/QWTic6HnLVN8EpnZsVv
2yQYnKc6F3A0gYGKUf+wzRR4XxG5Lx/TAC47Ijq0t82Axb3Cdre1OrjHfJV4
KDgDLc3f8dHROjCa3TYaLDYDEsmWo7e66gB++K91nZoBJkunJp+cOmDx99hn
enMGzD5tWXphVQdvmC0f2ubNgFb0qqBWby0kzMpxBV6dBT1UTXRU10DDj3bl
1sg5iD91VTa1rgryjQzlyxLnoLLYtG+osAriRwjJrMw5iKBr/ZROrwJH0vi2
N+Vz0GmMCkffVIGwNHlAfWwO2K5rLfzVrQKHgMNuZSfnYcrAllX+dyUIaXl/
yeqehwtbGqXeslfCw+8nxAJ46ZChpC2qYl8Onhtzz66K0IHd0j9S4045BMlk
UXcfosPREo++64z3cd57iZB0JTrIcIZExGPlsHRPgKXemlHq8TtRAVc5PN2y
OsBSQQdxSwOWLZ/K4JVmabTLg0UQ5cwwWxnDIbFZ+T+L2iXIfTZUY8yHgP99
0Hv71iX4dG8vCycLAj+LfpHnfUvg8Fb1Us10Kdj9eyIVN70Eu8J39j2oLwUl
yVzVXv5lsJ1S3uP+tBRo3iIuWlbLQH+uZOk7WgL8ius0OY4VmI9Kg5H6YvBP
/5qy6eIqrNw5mtHytRDEhM+ETV5dhW9t/ZvKPhYCHoy/6DBbBZUbv1VRSCHM
udSZJzitQrl906mRR4VgqP5TEBJWIX1QUnlMuhAODK2/dVxZhe3Upa74jG9Q
vee885+UNaBeda90zi4AcnC7egXrOgTMjxq2EfmwHj1mPkNZB8+llCf01nxY
SSG57tm5Dme3+6Xtq86H+dJDmU+k1mEr3+dtUen5MDTmte24/jpEeurQmV3y
oeGs7M/ElHX4Q7adSebNh/DZEGvvCxvAizRCyfpfgC81rSZRjYQl9Mcs7JnN
hW3/dqp/0CRhG1Jdy/nDubDl8qv62Csk7Ojdg4l61FxgWrZsjDQiYYNku+bU
2lyYVBNre2tLwtgCMTw9NheqR+JoruEkbOT07T182rngKBG6oDtCwuRnHXu1
c3OgM9ljH6sPGcvxeexPD8iGQW3OnvuvyZhNLNXdxzsbZhdD3nQHkbHMxdlE
Ebds4FZLW06LJmM1arWWVnez4fxIR/OVXDJm6DDPr6WSDV8lJF3jaGSsP/jI
WYu1LAhP7Wg/LcuEJWLFFokPsuB62kFvu19MWMkQ67lNrZmwc8/6D/dhJiz/
o3n5vYpM6A1q2fV6ggnTtBc/1JuXCYYuj7OT6UwY747pqPawTDBQr+7r42DG
3PtebgLjTNAdNpNRlWHGLi8HFv+ayICL+yL6d3kxYxuZjkcOb82AM+83na4S
ZcHu7uvu22efBr8DvtH3e7FiCpWsJu4mKZC3z31cpI8Ns7Na6BT8kQT9CU7H
WHazY+fcm4XUKxKgVULFY+7xZozbZSC7rSse9kSqCNeWcGDKeouxkepxEPrp
2vLJNU6MuPDL8/WfaGhDZ3sK9LiwoeIAFW/993CvzJhzdzg3pmRMbf55JgLM
Y3tCi6O5sYuFLQrMxyPA0E13r1ECNzZQo6J4RDIC1E5qnHyfyY11lGmNvROI
gP3ZsncEqrixVioRHkoPh764LRX/zXJjKYpiKbbZ4XDRo9CF+yIFS7pEab17
MBzOGcuyZV6mYEbdbmVCe8Ph1JmsQE19CpbAIf2nUyAcJJaSUvzNKZj4pFyp
zuZwYLEN7ORwpmAGHMes3EfCoMTE8timjxTs9PnLClHJYZCn8LvkYxoF83w0
IBYUGwbpu4zVz2dTMC8rXgH/sDCI7NIxfVFMwQaj2rb4vwwDR23GSGyjYL77
O4aL7oSBzeFCgYRuClbkkd/QYhIGtzhlk87SKNiZCoWkcf0wuFJzsPjZKMOe
F0NikmphoJ6UpCoyRcEWngzUX8LCAHsu3IbPUzDDs7IGLvJhcMLsvbHpCgXr
L9JpSJYOg//7/8P+///vfwBc30qm
     "]]}, 
   {RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], Opacity[
    1.], LineBox[CompressedData["
1:eJw1mHk01V3Ux69Zwu2RkkgoQ6inkjHal2RWKpUiiUxFpvCghJApitJEyJCp
i0wpnJ+Qeegar5Kp4ZruL9MN1/Ded633/eusz1r7nHX2d5/vWXttKTv3Mw7s
BAJBjI1A+N/1AGyf3NggYnIhf/QkJMhgcV0ucnmFiPXvNDnEI0eGgCfqsvOL
RCypSmtt/l8yNE5dsvs1RcS0UKEWVYcMts9SB9v6iZhqUelTigMZHs3Ktz4v
JGJvdQzehZLJsJBxtED5ChE7/fm6K023ED7wXr3pVCOIxTKIt+S0iuBue99U
7VUBLFtc0VQhvxgmskRF1MT5sZM9CXVmB0pAYdhwE3GIDxt0feH5ZrAUMrdm
LFBDN2E3Mic6sqfKoc6q5QECXmxCw3i8V7ISxFPFJZ+NcmMfKoUvpo19BJPj
Vv7dUVyYXm5s3tFNCHYzHYP3GXFi4V49T4w1a2E/OSXYa4YdE4mQKn/mWAcj
t9IP9OWwYRN2rk7prQ0QpOJjOJTBhj19VIS6extAnGFo9+MVG1ZU7VrHP9IA
ln5/nsw9ZsP8RtLHEhYa4Is/rAmGsGFG9y23/9n1GT4FDbUaXGTD5m5c2jbv
+RkyI0VdKjexYfU2bvd0dzeCS/KjjOTrBOy/n5xa5g+bYIHXOErNgYCZb4mY
sU1pgmAfdnfKFQL2LqBMIzCvCZJOemvyniNgMrc46tvqm+AT4dwXbyBgB/de
iP+x3AQ7HXZsGAsTsDJZtkS6QzM8sTpWNhe4gSY7pmLyTrSAZdDFPQXS60jy
182pTpk22Iyp17/YsY4iesWEQ460QQ3bDocowXXUEi+cp3W8DWTC+7IdV9bQ
E/tbpG7bNpiLOasgRVlDRpoGfHeT2yD2udmhpJA1tF3buX5gWztgpToQPLKK
qjtXh7q3dID81D4ri1dMJBW2LXxVuAvKwiiJ1xOZ6LipkkmYRBfoSgS2BUcy
UY6I5LOt8l1gfbpN+60XE70qOjNqeLQLHlbclOQ2YCLzK3PTgnZdsHyv5GcF
voK07R4NChZ3QYu4toeYzgoqPNeTK3juC9w4aR4+Nr6E6iSv3c8op8C68cEe
ReoSquw6dULwEwUSDLbs8elYQlf5f3OFtlPgA6kL46lcQpx/VgvCflCAV9l8
VSluCfWWj1/7K9QN2SLm3n7qS+hGvZDoc89uGBs5dZX/wV9k3Jl9lEO1By55
n9JWVWUgPzmDon1DvRAzkHJVXYmBgvgiG6Qme6FKezpcU5qBmB1+Dnv+9oIE
b1THMUEGat2TfuX4P30wmvzpisGvRSQt1bHAONEHzp9VQi2TFlEuYfDK4eI+
8BUVbwxgLKCF1Ndc8vH9kIho5qh8HmGPlVP9rlPh1W2CwET+PBqtaW84H0iF
HI0dzULp82jZRIP7eCwVakoMdBxj5hEu8HROl0wF2pvsQwK28+itKN1odI4K
xx5eE7LaNI8C+5+ffxU0CBO2I90Mqzlke9G1wDjtKwD7wPkDbLNIT+FkbK3A
d7g+8Mkja34acUk9Wsl7Ngpaz0N2fPjvN2LynKDYdfwA56HZrxdFxpCgQFp6
P9CgxP0blz7HGCrbubVixJQG62yN/x7GR5GS3w/SxEUaPJZNvsfXOIq8VbnM
l7xpgLnrK370HUU7HLmyN+fQYAf7ywDxvhEkdmHp6x/iBDw0izxcs28YSTpX
Vi4NTcCpEaW3TXsHkb+p1alg1ykwvSZnC8tUlPqw2S7GfwqMaVJby9upaLNG
1kpsxBTo/dnun+FLRdpv+Lm806ZAncB2IqhpAH1aPPOioXsKpCV7vym79qPq
swLcl9WnYfHKHf7Ukh6kbCEd8M/SNMyN+6FtkT0oQf8fl8ccM/DHycsr1roH
hfpmb+MlzsCku9OAP1cPUgirnMJkZuD73dOZFhe60V6it7P3mRloSpXR4mN+
QYNqNPK7nBlQRS/7J8M6UDRhz84LBnTYyc858Uq5A301XbFhnqLDmqXrypmx
dtR/2Mn6gSUd6ue0dn2EdjRsZbHg70KHs7Lf7WJWWlGy1EizZzQdPB5I0hU9
mpFSSIBWTRMdLAYjN4YlmtHNN40L+7vooCY3u+VxexMKNDlkHNVPh3WsVnlV
oQmVM93N2H/RIXbBzr/t52dkKdpvssKGQ65VFqebVT1697br028VHGJzBLZL
8dUjRRmddXYtHDwWfeR639chwzMpTB5dHNTj9Y21t9ehboHVzg4zHD5/+h0v
8KUWhex0WDO9hkMe0Ty99m4tokWxXStxweGB9ft3PgdqUaeQcjW7Ow4WjMje
oRgMPeS7+svKn7X/rIWq9YMapPBxLWQxlhU/s/di/WQ1clHrpaU9YsXfXwxU
MqxGnO1i0UeScBj/kFS7yl6FLBPU5YRTcSBIU01T/qtEm1fWHr4rxKHxY647
V997tF12+EFKCQ5x5wIS3JTfo82KJFmPChzEo8QGtOnlyM1serS1BgeNP5ft
v9uXIt+TvW+s2nBw/ZA54VdfgthDNfm8OnF4FTbl/o9MCeqT8rjtQcGBQ9T/
tt6vYlR4ZeGizAAOquOI/bt+MTL2ksz/OoiD81vuKL83RegWB7XffwiH5/fk
Fc99ISOeOwaP3cZwWBcpV4tMeIt+/KcuUPcDB7uC43ofzxagOIn9Rwi/cVDq
tbks1ZeLlK6a2B+awuGRy7SLxdMcdGhhl6XMDA6MNX/f+5Zv0Fmmfg8BxwGT
fRI/M5iJXtSVkW7O4SDzUTpZMjkD2YOkFGEBh+hTRTlnL79GEVaKLv6LOODj
2mUREumIS5hDaJCBw5s/Sfe8USqaSNHbJr3Eqrfou2SqTgoSK153Pb2Mg6Zu
exnUv0DK7cVbnVdw4LxB68jSf4aU/c05HZk4dCRy0DY3P0FbgkDTZBWHZ1US
bF4miUhO/HyxyBorv58aOwfaH6LAo98dWll8QTzP8XZaLOpJemzhuI7DpNvI
M8Os+yh1NCeAxuIgtL1VOC8UpdGqe09v4CAt/WKl9nIguv52wv41i40auK9c
cvJE2bfL9gyyGEtvF2itu4wEN3YJL7PYMGMDy1bUA50TiUc2WBz22vH3+7/O
AB1LYdMsFiPrE/W8fKE6MourlsUBXG69544Eg1qnfcVtFv8qUc2I+jccimw8
X+5m8Vk7gme1QjScMrxbksu6H7al5disTDy8GypmE2Nx0rBRRMKWBGB/lhbp
zcqvrfyAm1LcYxhdPmdSwtKDPW6rxefNT2GT6ymzQZZe6g5LmrZRz6Fd+E4C
jaWnu9aQ1Ap3MlzL09n9naV39tZPvI/DXgEbyWvpA6seN326/24PSYMOgfN7
7v7FIdP0/fei9deg4/cgvJpVX3ivQLmjnQmzzDfNavM4DO5JaTC+nQX8f9n9
n8/isIUZUvBj+Q3MxY+pbaXjkO+4kFqsngu7V0ouKE7joE9xTAzyy4P7FTYc
ipM43Mk1CdixWABWTzP3jP7EYce2GrefymSwz+cueDqOQ2nwwavvvAqhWajz
t8ooDov5fHtll4qA9mQga983HBTv4XyYTwl0PfCKwrtwsD66lfetdCl4dse0
p7ez/pd5Va4XnaWw8noXrtaCA90+aMN7Xzm4/00Q3lKHwzs9gQW5offwIDlU
wIvl3/HVg7PC0ZXgpiJuc4Llb+EyCzqb2gfwPOF9YC0PBx+ZZNrXhx+BsLcv
d+drlr7cikPxejUwEH4ylxiPg0vNycHbszUgWTj5aFM0Di98vfpdXiHoP9xY
Qwtjvbf5bzSp2xgwtOVOXwjAYdh+7ziztxaOPOB0WrLH4bVeSX9hVD2YyiVb
XDuMwwfJE+dlR+rBKC5JElPEoXu1rydZtQF8JI23c8iw/FG28iVqvAHSsq7x
gAgODjK6bfbajWChxEtJWKGDPPcXTGS2GVS0bdubq+hAGrODOIMW6AwYZuaU
0uFizUI156sW0G0rT/MsoEO0746Ps8atcNVXYRq9pMP0rytlrZltMJGpM2vh
T4fCRnruXctO6HWtaw47SIfvNgJ3/s3thFv9zxhFcnQQZCiaDy+z2MKtrlWC
DqoJRwJTzbrgt5FTSzs/HR5jnlIbdFbfdeCaFw9tBk7vmnarOUSBPzwLJR4v
Z6Clb4Rbu6IHngzqLor9nYZqoxZ1zU9USDedz1iOmgK6ypG6+Rkq3JXnlJ++
PQUS0q/M3ooOQv6uFmKv+xQErXjaSXoOwtj95Jroc1NAyheN5Zb8ChX1J+9E
S05BPb/LMOXON5AlxWxmL5uElk6eiBsaw/CigxrBHJiAqbzXLictx+CiYEJT
zjYapLw98HKh7DdUHU9hC8r5AVVXI7YdrJiGfF3vH4OVoxAg/Lz3lNAsqOyN
cyqbGYIbHOeHPcVnIfT9R/KXwSGwnhOaSJSdBQfGsfSVxiGArtjVfs1ZWHqX
f9vr9RBwxgbvsbWfBUbh5YrlC0MQx+Hs6VE6CyPIzsy+4RtkzKkKPDo/B0rD
sz33c79CW1ePHuXFPFxBief5k6hQbn1JvTZzHjb99v8WEsHSiTasWESeB5rT
dWd2Pyr4EqaE4j7NQ5y/UuBeSypIHmQbNZqch/rZI7SpnVTwid9/p1ZjAYj7
VMmhrwdA4mREWdHAAoQeHLbNquwHz1Y12fjtDCiapx2z5e6D0I35kPNSDDCa
9iowXeqFBOWioV1KDIDwDS09Vn9a+lL+SYEOA7KUDmZbdPTC0g1RzhZXBozr
+rxbTuqFIH7mKGcdA5xBPsRMoReizWpSAtz/goPLn7djlj2Q2am3zbFpCWD0
+fpGNwV2vEx4eYuyBFE8unOazRSIcRyRuvdtCdzH7w3eqaGA11rggbQ/S/AP
Q8ZDLJcCOoolBoM7loE53cJRFkSB7xFSASddliE2Rt/XWpF13rH176p8K6Dl
Zxt6JvoLxBZU5PCYMmGXPtMx6XIXyEpqPaWfZ8LKXKn38NkuwBKx8N6rTOAw
C5c/aNwF8wHNdhl+TBAy8HOnqXbBJaOv4pDBhBhriYqSLV2w79f6Q98VJtRk
19bbxHTC5936/j9zVoH3WE2N2PUOYEvsMarjWocCBWlucfM2WE+ZtJslroMT
7dWx/TptsJJDuL175zpUeXTLGR1ug4UaJXLggXVQ7sZoacJt8GsyTOjIhXXQ
buxW6qS2QpuuytfMnHXQdNoqstuxFZ7NPXGNMNmAKQ1CpXlEC4jk5jdmGhJI
jblyLkfGm6DvTfAerkg2krDszgsuxM9gla8Q4TXGTqpfD6Q2+teB1kueow0y
nKSz7GNV2M1aGI+vZOwN4yL9peoIZIogKN1zd0rqGzfJ4G5ZVTD9I4xk+B3m
3MVLSk++/0jTrhIo8ieC5//bRBLbtj9mlFkOu1+ckGyq5iPpmfsIRfKUQVLe
xWWN1c2k8sJEZ/r5EuhGutT35wRIMYli2fnvi+FGrc3mXc8ESVvimBRDuyKw
S6UmVaUIkpSbbUTfWBfBpTsW0tYZgqQ49xoV3gtFYKhhrPGSLEj6ZMUzQTEp
gr3FKk6iDYKkjVJJ7miVIviWxl+3bU6QZLAOoXc2FYFp8IcAQVMi6ato9fZB
x0I4bqPCTT5NJOWMlcqfsC4ETa2iR2YXiKStmy75lJ0uBPmlrJxYOyKpD2kW
Z2kVAqfHoz4+fyJpnLv9XL1QIVRfcT7Mk00kCRp43fyLyFCqPV6dnU8kWdwz
TyktI0OBmI2RfjGRVFqlQ/HNJ8OL/rO24VVEkm5AvA7fUzL4nmJZvptIWtSt
++Z6kwxu+z+IZgwQSSTC1vHT18hwbbNKlu53Iml30s6+o5fIcKZRoSpkgkgS
Vqb4SOiTwSgry0AKJ5JaKwzFRbTIQLon2Y0tEEls/JwZwofJoHb1pY3tCpGU
dnOGW0SeDP83nyH9/3zmfwD6ZhVU
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  Method->{
   "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None},
  PlotRange->{{-4, 4}, {0., 0.9999996210919236}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.8227054913300934`*^9, {3.8227055505187144`*^9, 3.8227055848419333`*^9}, 
   3.8227108492393403`*^9, 3.8227273657580204`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Series", "[", 
  RowBox[{
   RowBox[{"tf", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "1"}], "/", 
      RowBox[{"(", 
       RowBox[{"1", "+", 
        RowBox[{"x", "^", "2"}]}], ")"}]}], ",", "2"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "Infinity", ",", "4"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8227056040720835`*^9, 3.8227056382251945`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
   FractionBox["1", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "+", 
       SuperscriptBox["\[ExponentialE]", "2"]}], ")"}], " ", 
     SuperscriptBox["x", "2"]}]], "+", 
   FractionBox[
    RowBox[{
     RowBox[{"-", "2"}], "-", 
     SuperscriptBox["\[ExponentialE]", "2"]}], 
    RowBox[{"2", " ", 
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{"1", "+", 
        SuperscriptBox["\[ExponentialE]", "2"]}], ")"}], "2"], " ", 
     SuperscriptBox["x", "4"]}]], "+", 
   InterpretationBox[
    SuperscriptBox[
     RowBox[{"O", "[", 
      FractionBox["1", "x"], "]"}], "5"],
    SeriesData[Parallel`Preferences`Private`x, 
     DirectedInfinity[1], {}, 2, 5, 1],
    Editable->False]}],
  SeriesData[Parallel`Preferences`Private`x, 
   DirectedInfinity[
   1], {(1 + E^2)^(-1), 0, Rational[1, 2] (-2 - E^2) (1 + E^2)^(-2)}, 2, 5, 1],
  Editable->False]], "Output",
 CellChangeTimes->{{3.822705613607503*^9, 3.8227056389852157`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"LaplaceTransform", "[", 
  RowBox[{
   RowBox[{"1", "/", 
    RowBox[{"(", 
     RowBox[{"1", "+", 
      RowBox[{"x", "^", "2"}]}], ")"}]}], ",", "x", ",", "w"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.822708411914157*^9, 3.8227084247879868`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"CosIntegral", "[", "w", "]"}], " ", 
   RowBox[{"Sin", "[", "w", "]"}]}], "+", 
  RowBox[{
   FractionBox["1", "2"], " ", 
   RowBox[{"Cos", "[", "w", "]"}], " ", 
   RowBox[{"(", 
    RowBox[{"\[Pi]", "-", 
     RowBox[{"2", " ", 
      RowBox[{"SinIntegral", "[", "w", "]"}]}]}], ")"}]}]}]], "Output",
 CellChangeTimes->{3.8227084250562773`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Integrate", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"1", "/", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"k", "-", "p"}], ")"}], "^", "2"}], "+", 
       RowBox[{"c", "^", "2"}]}], ")"}]}], "*", 
    RowBox[{"1", "/", 
     RowBox[{"p", "^", "2"}]}]}], ",", 
   RowBox[{"{", 
    RowBox[{"p", ",", "pf", ",", "Infinity"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8227094036780787`*^9, 3.822709426998632*^9}}],

Cell[BoxData["$Aborted"], "Output",
 CellChangeTimes->{3.822709435637659*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Integrate", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"1", "/", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"k", "+", "p"}], ")"}], "^", "2"}], "+", 
       RowBox[{"1", "^", "2"}]}], ")"}]}], "*", 
    RowBox[{"1", "/", 
     RowBox[{"p", "^", "2"}]}]}], ",", "p"}], "]"}]], "Input",
 CellChangeTimes->{{3.8227094036780787`*^9, 3.8227094480391665`*^9}, {
  3.82270965654217*^9, 3.822709656664837*^9}, {3.8227096956264677`*^9, 
  3.8227097101108627`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{"-", "1"}], "-", 
   SuperscriptBox["k", "2"], "+", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      SuperscriptBox["k", "2"]}], ")"}], " ", "p", " ", 
    RowBox[{"ArcTan", "[", 
     RowBox[{"k", "+", "p"}], "]"}]}], "-", 
   RowBox[{"2", " ", "k", " ", "p", " ", 
    RowBox[{"Log", "[", "p", "]"}]}], "+", 
   RowBox[{"k", " ", "p", " ", 
    RowBox[{"Log", "[", 
     RowBox[{"1", "+", 
      SuperscriptBox["k", "2"], "+", 
      RowBox[{"2", " ", "k", " ", "p"}], "+", 
      SuperscriptBox["p", "2"]}], "]"}]}]}], 
  RowBox[{
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"1", "+", 
      SuperscriptBox["k", "2"]}], ")"}], "2"], " ", "p"}]]], "Output",
 CellChangeTimes->{
  3.822709448734262*^9, 3.822709658968093*^9, {3.822709697072302*^9, 
   3.8227097104370413`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FullSimplify", "[", "%", "]"}]], "Input",
 CellChangeTimes->{{3.8227095890466638`*^9, 3.8227095936439915`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{"c", "-", "k"}], ")"}], " ", 
    RowBox[{"(", 
     RowBox[{"c", "+", "k"}], ")"}], " ", "p", " ", 
    RowBox[{"ArcTan", "[", 
     FractionBox[
      RowBox[{"k", "-", "p"}], "c"], "]"}]}], "-", 
   RowBox[{"c", " ", 
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox["c", "2"], "+", 
      SuperscriptBox["k", "2"], "+", 
      RowBox[{"k", " ", "p", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"Log", "[", 
          RowBox[{
           SuperscriptBox["c", "2"], "+", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"k", "-", "p"}], ")"}], "2"]}], "]"}], "-", 
         RowBox[{"2", " ", 
          RowBox[{"Log", "[", 
           RowBox[{"-", "p"}], "]"}]}]}], ")"}]}]}], ")"}]}]}], 
  RowBox[{"c", " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox["c", "2"], "+", 
      SuperscriptBox["k", "2"]}], ")"}], "2"], " ", "p"}]]], "Output",
 CellChangeTimes->{3.822709593912927*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"an", "[", 
   RowBox[{"x_", ",", "n_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{"1", "/", "2"}], "/", "\[Pi]"}], "*", 
   RowBox[{"n", "/", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"n", "/", "2"}], ")"}], "^", "2"}], "+", 
      RowBox[{"x", "^", "2"}]}], ")"}]}]}]}]], "Input",
 CellChangeTimes->{{3.822710584590319*^9, 3.822710615126859*^9}, {
  3.8227114035590067`*^9, 3.822711406102501*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"an", "[", 
  RowBox[{"x", ",", "n"}], "]"}]], "Input",
 CellChangeTimes->{{3.82271061631169*^9, 3.8227106184963703`*^9}}],

Cell[BoxData[
 FractionBox["n", 
  RowBox[{"2", " ", "\[Pi]", " ", 
   RowBox[{"(", 
    RowBox[{
     FractionBox[
      SuperscriptBox["n", "2"], "4"], "+", 
     SuperscriptBox["x", "2"]}], ")"}]}]]], "Output",
 CellChangeTimes->{3.822710618843954*^9, 3.822711408378039*^9, 
  3.822727388136547*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Integrate", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"an", "[", 
       RowBox[{
        RowBox[{"k", "+", "p"}], ",", "n"}], "]"}], "/", 
      RowBox[{"p", "^", "2"}]}], ",", "p"}], "]"}], "*", "2", "*", 
   RowBox[{"\[Pi]", "/", "4"}]}], "//", "Simplify"}]], "Input",
 CellChangeTimes->{{3.822710620694624*^9, 3.8227106434788437`*^9}, {
  3.822711414669593*^9, 3.8227114504729023`*^9}, {3.822711516070405*^9, 
  3.822711516206392*^9}}],

Cell[BoxData[
 RowBox[{"-", 
  RowBox[{
   FractionBox["1", 
    RowBox[{
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{
        RowBox[{"4", " ", 
         SuperscriptBox["k", "2"]}], "+", 
        SuperscriptBox["n", "2"]}], ")"}], "2"], " ", "p"}]], 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"2", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"-", "4"}], " ", 
         SuperscriptBox["k", "2"]}], "+", 
        SuperscriptBox["n", "2"]}], ")"}], " ", "p", " ", 
      RowBox[{"ArcTan", "[", 
       FractionBox[
        RowBox[{"2", " ", 
         RowBox[{"(", 
          RowBox[{"k", "+", "p"}], ")"}]}], "n"], "]"}]}], "+", 
     RowBox[{"n", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"4", " ", 
         SuperscriptBox["k", "2"]}], "+", 
        SuperscriptBox["n", "2"], "+", 
        RowBox[{"8", " ", "k", " ", "p", " ", 
         RowBox[{"Log", "[", "p", "]"}]}], "-", 
        RowBox[{"4", " ", "k", " ", "p", " ", 
         RowBox[{"Log", "[", 
          RowBox[{
           RowBox[{"4", " ", 
            SuperscriptBox["k", "2"]}], "+", 
           SuperscriptBox["n", "2"], "+", 
           RowBox[{"8", " ", "k", " ", "p"}], "+", 
           RowBox[{"4", " ", 
            SuperscriptBox["p", "2"]}]}], "]"}]}]}], ")"}]}]}], 
    ")"}]}]}]], "Output",
 CellChangeTimes->{{3.822710633092724*^9, 3.822710644650977*^9}, {
   3.8227114106209855`*^9, 3.822711420249496*^9}, 3.822711450830872*^9, 
   3.8227115167523246`*^9, 3.822727390027131*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Simplify", "[", "%", "]"}]], "Input",
 CellChangeTimes->{{3.8227106477945347`*^9, 3.8227106512070246`*^9}}],

Cell[BoxData[
 RowBox[{"-", 
  RowBox[{
   FractionBox["1", 
    RowBox[{
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{
        RowBox[{"4", " ", 
         SuperscriptBox["k", "2"]}], "+", 
        SuperscriptBox["n", "2"]}], ")"}], "2"], " ", "p"}]], 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"2", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"-", "4"}], " ", 
         SuperscriptBox["k", "2"]}], "+", 
        SuperscriptBox["n", "2"]}], ")"}], " ", "p", " ", 
      RowBox[{"ArcTan", "[", 
       FractionBox[
        RowBox[{"2", " ", 
         RowBox[{"(", 
          RowBox[{"k", "+", "p"}], ")"}]}], "n"], "]"}]}], "+", 
     RowBox[{"n", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"4", " ", 
         SuperscriptBox["k", "2"]}], "+", 
        SuperscriptBox["n", "2"], "+", 
        RowBox[{"8", " ", "k", " ", "p", " ", 
         RowBox[{"Log", "[", "p", "]"}]}], "-", 
        RowBox[{"4", " ", "k", " ", "p", " ", 
         RowBox[{"Log", "[", 
          RowBox[{
           RowBox[{"4", " ", 
            SuperscriptBox["k", "2"]}], "+", 
           SuperscriptBox["n", "2"], "+", 
           RowBox[{"8", " ", "k", " ", "p"}], "+", 
           RowBox[{"4", " ", 
            SuperscriptBox["p", "2"]}]}], "]"}]}]}], ")"}]}]}], 
    ")"}]}]}]], "Output",
 CellChangeTimes->{
  3.822710651459313*^9, {3.822711422993836*^9, 3.822711437398667*^9}, 
   3.822727395671734*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"f", "=", 
  RowBox[{
   RowBox[{"2", "*", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", "*", 
       RowBox[{"k", "^", "2"}]}], "-", 
      RowBox[{"n", "^", "2"}]}], ")"}], "*", "p", "*", 
    RowBox[{"ArcTan", "[", 
     RowBox[{
      RowBox[{"2", "/", "n"}], "*", 
      RowBox[{"(", 
       RowBox[{"k", "+", "p"}], ")"}]}], "]"}]}], "-", 
   RowBox[{"n", "*", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", "*", 
       RowBox[{"k", "^", "2"}]}], "+", 
      RowBox[{"n", "^", "2"}]}], ")"}]}], "-", 
   RowBox[{"4", "*", "n", "*", "k", "*", "p", "*", 
    RowBox[{"Log", "[", 
     RowBox[{
      RowBox[{"p", "^", "2"}], "/", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"n", "^", "2"}], "+", 
        RowBox[{"4", "*", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"k", "+", "p"}], ")"}], "^", "2"}]}]}], ")"}]}], 
     "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.8227117400036354`*^9, 3.822711792250403*^9}, {
  3.8227119289443827`*^9, 3.8227119408476353`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"-", "n"}], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"4", " ", 
      SuperscriptBox["k", "2"]}], "+", 
     SuperscriptBox["n", "2"]}], ")"}]}], "+", 
  RowBox[{"2", " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"4", " ", 
      SuperscriptBox["k", "2"]}], "-", 
     SuperscriptBox["n", "2"]}], ")"}], " ", "p", " ", 
   RowBox[{"ArcTan", "[", 
    FractionBox[
     RowBox[{"2", " ", 
      RowBox[{"(", 
       RowBox[{"k", "+", "p"}], ")"}]}], "n"], "]"}]}], "-", 
  RowBox[{"4", " ", "k", " ", "n", " ", "p", " ", 
   RowBox[{"Log", "[", 
    FractionBox[
     SuperscriptBox["p", "2"], 
     RowBox[{
      SuperscriptBox["n", "2"], "+", 
      RowBox[{"4", " ", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"k", "+", "p"}], ")"}], "2"]}]}]], "]"}]}]}]], "Output",
 CellChangeTimes->{3.8227117929146566`*^9, 3.8227119422183914`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"f", "=", 
  RowBox[{
   RowBox[{"f", "/", "p"}], "/", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", "*", 
       RowBox[{"k", "^", "2"}]}], "+", 
      RowBox[{"n", "^", "2"}]}], ")"}], "^", "2"}]}]}]], "Input",
 CellChangeTimes->{{3.822711795504298*^9, 3.8227118036579247`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{
    RowBox[{"-", "n"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "+", 
      SuperscriptBox["n", "2"]}], ")"}]}], "+", 
   RowBox[{"2", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "-", 
      SuperscriptBox["n", "2"]}], ")"}], " ", "p", " ", 
    RowBox[{"ArcTan", "[", 
     FractionBox[
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{"k", "+", "p"}], ")"}]}], "n"], "]"}]}], "-", 
   RowBox[{"4", " ", "k", " ", "n", " ", "p", " ", 
    RowBox[{"Log", "[", 
     FractionBox[
      SuperscriptBox["p", "2"], 
      RowBox[{
       SuperscriptBox["n", "2"], "+", 
       RowBox[{"4", " ", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{"k", "+", "p"}], ")"}], "2"]}]}]], "]"}]}]}], 
  RowBox[{
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "+", 
      SuperscriptBox["n", "2"]}], ")"}], "2"], " ", "p"}]]], "Output",
 CellChangeTimes->{3.822711804054032*^9, 3.822711945646871*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"D", "[", 
   RowBox[{"f", ",", "p"}], "]"}], "//", "FullSimplify"}]], "Input",
 CellChangeTimes->{{3.8227118066891117`*^9, 3.822711815407117*^9}}],

Cell[BoxData[
 FractionBox["n", 
  RowBox[{
   SuperscriptBox["p", "2"], " ", 
   RowBox[{"(", 
    RowBox[{
     SuperscriptBox["n", "2"], "+", 
     RowBox[{"4", " ", 
      SuperscriptBox[
       RowBox[{"(", 
        RowBox[{"k", "+", "p"}], ")"}], "2"]}]}], ")"}]}]]], "Output",
 CellChangeTimes->{3.8227118159327574`*^9, 3.8227118793012395`*^9, 
  3.8227119477497406`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"%", "-", 
   RowBox[{"2", "*", 
    RowBox[{"\[Pi]", "/", "4"}], "*", 
    RowBox[{
     RowBox[{"an", "[", 
      RowBox[{
       RowBox[{"k", "+", "p"}], ",", "n"}], "]"}], "/", 
     RowBox[{"p", "^", "2"}]}]}]}], "//", "FullSimplify"}]], "Input",
 CellChangeTimes->{{3.8227118178911524`*^9, 3.822711871378945*^9}}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{{3.822711844031114*^9, 3.8227118817659903`*^9}, 
   3.8227119504711027`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"fint", "[", 
   RowBox[{"k_", ",", "kc_", ",", "n_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"n", "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"4", "*", 
            RowBox[{"k", "^", "2"}]}], "+", 
           RowBox[{"n", "^", "2"}]}], ")"}]}], "-", 
        RowBox[{"2", "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"4", "*", 
            RowBox[{"k", "^", "2"}]}], "-", 
           RowBox[{"n", "^", "2"}]}], ")"}], "*", "kc", "*", 
         RowBox[{"ArcTan", "[", 
          RowBox[{
           RowBox[{"2", "/", "n"}], "*", 
           RowBox[{"(", 
            RowBox[{"k", "+", "kc"}], ")"}]}], "]"}]}], "+", 
        RowBox[{"4", "*", "n", "*", "k", "*", "kc", "*", 
         RowBox[{"Log", "[", 
          RowBox[{
           RowBox[{"kc", "^", "2"}], "/", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"n", "^", "2"}], "+", 
             RowBox[{"4", "*", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{"k", "+", "kc"}], ")"}], "^", "2"}]}]}], ")"}]}], 
          "]"}]}]}], ")"}], "/", "kc"}], "/", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"n", "^", "2"}], "+", 
        RowBox[{"4", "*", 
         RowBox[{"k", "^", "2"}]}]}], ")"}], "^", "2"}]}], "*", 
    RowBox[{"2", "/", "\[Pi]"}]}], "+", 
   RowBox[{"2", "*", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"4", "*", 
        RowBox[{"k", "^", "2"}]}], "-", 
       RowBox[{"n", "^", "2"}]}], ")"}], "/", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"4", "*", 
         RowBox[{"k", "^", "2"}]}], "+", 
        RowBox[{"n", "^", "2"}]}], ")"}], "^", "2"}]}]}]}]}]], "Input",
 CellChangeTimes->{{3.8227123336275983`*^9, 3.822712436540536*^9}, {
  3.822712681082164*^9, 3.8227126974832025`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"fint", "[", 
   RowBox[{
    RowBox[{"-", "1"}], ",", "2", ",", "4"}], "]"}], "//", "N"}]], "Input",
 CellChangeTimes->{{3.8227124383894606`*^9, 3.822712444250539*^9}, {
  3.8227132206984987`*^9, 3.822713243104162*^9}}],

Cell[BoxData["0.06235601125295223`"], "Output",
 CellChangeTimes->{
  3.822712445385111*^9, {3.8227127087573338`*^9, 3.8227127393139095`*^9}, {
   3.822713204674821*^9, 3.8227132434741893`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData["f"], "Input",
 CellChangeTimes->{3.8227129589845123`*^9}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{
    RowBox[{"-", "n"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "+", 
      SuperscriptBox["n", "2"]}], ")"}]}], "+", 
   RowBox[{"2", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "-", 
      SuperscriptBox["n", "2"]}], ")"}], " ", "p", " ", 
    RowBox[{"ArcTan", "[", 
     FractionBox[
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{"k", "+", "p"}], ")"}]}], "n"], "]"}]}], "-", 
   RowBox[{"4", " ", "k", " ", "n", " ", "p", " ", 
    RowBox[{"Log", "[", 
     FractionBox[
      SuperscriptBox["p", "2"], 
      RowBox[{
       SuperscriptBox["n", "2"], "+", 
       RowBox[{"4", " ", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{"k", "+", "p"}], ")"}], "2"]}]}]], "]"}]}]}], 
  RowBox[{
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "+", 
      SuperscriptBox["n", "2"]}], ")"}], "2"], " ", "p"}]]], "Output",
 CellChangeTimes->{3.8227129597440357`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"f", "/.", 
  RowBox[{"p", "\[Rule]", 
   SubscriptBox["k", "c"]}]}]], "Input",
 CellChangeTimes->{{3.822713020498342*^9, 3.8227130246575003`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{
    RowBox[{"-", "n"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "+", 
      SuperscriptBox["n", "2"]}], ")"}]}], "+", 
   RowBox[{"2", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "-", 
      SuperscriptBox["n", "2"]}], ")"}], " ", 
    RowBox[{"ArcTan", "[", 
     FractionBox[
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{"k", "+", 
         SubscriptBox["k", "c"]}], ")"}]}], "n"], "]"}], " ", 
    SubscriptBox["k", "c"]}], "-", 
   RowBox[{"4", " ", "k", " ", "n", " ", 
    RowBox[{"Log", "[", 
     FractionBox[
      SubsuperscriptBox["k", "c", "2"], 
      RowBox[{
       SuperscriptBox["n", "2"], "+", 
       RowBox[{"4", " ", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{"k", "+", 
           SubscriptBox["k", "c"]}], ")"}], "2"]}]}]], "]"}], " ", 
    SubscriptBox["k", "c"]}]}], 
  RowBox[{
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "+", 
      SuperscriptBox["n", "2"]}], ")"}], "2"], " ", 
   SubscriptBox["k", "c"]}]]], "Output",
 CellChangeTimes->{3.822713025285902*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"fint", "[", 
  RowBox[{"k", ",", 
   SubscriptBox["k", "c"], ",", "n"}], "]"}]], "Input",
 CellChangeTimes->{{3.822712484875309*^9, 3.8227125177874603`*^9}}],

Cell[BoxData[
 RowBox[{
  FractionBox[
   RowBox[{"2", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "-", 
      SuperscriptBox["n", "2"]}], ")"}]}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "+", 
      SuperscriptBox["n", "2"]}], ")"}], "2"]], "+", 
  FractionBox[
   RowBox[{"2", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"n", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"4", " ", 
          SuperscriptBox["k", "2"]}], "+", 
         SuperscriptBox["n", "2"]}], ")"}]}], "-", 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"4", " ", 
          SuperscriptBox["k", "2"]}], "-", 
         SuperscriptBox["n", "2"]}], ")"}], " ", 
       RowBox[{"ArcTan", "[", 
        FractionBox[
         RowBox[{"2", " ", 
          RowBox[{"(", 
           RowBox[{"k", "+", 
            SubscriptBox["k", "c"]}], ")"}]}], "n"], "]"}], " ", 
       SubscriptBox["k", "c"]}], "+", 
      RowBox[{"4", " ", "k", " ", "n", " ", 
       RowBox[{"Log", "[", 
        FractionBox[
         SubsuperscriptBox["k", "c", "2"], 
         RowBox[{
          SuperscriptBox["n", "2"], "+", 
          RowBox[{"4", " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"k", "+", 
              SubscriptBox["k", "c"]}], ")"}], "2"]}]}]], "]"}], " ", 
       SubscriptBox["k", "c"]}]}], ")"}]}], 
   RowBox[{
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       RowBox[{"4", " ", 
        SuperscriptBox["k", "2"]}], "+", 
       SuperscriptBox["n", "2"]}], ")"}], "2"], " ", "\[Pi]", " ", 
    SubscriptBox["k", "c"]}]]}]], "Output",
 CellChangeTimes->{
  3.822712487989609*^9, 3.822712518776237*^9, {3.8227127015409613`*^9, 
   3.822712724754073*^9}, 3.822712924222316*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"N", "[", "%", "]"}]], "Input",
 CellChangeTimes->{{3.8227124462529054`*^9, 3.822712447350293*^9}}],

Cell[BoxData[
 RowBox[{"-", "0.10916454869905547`"}]], "Output",
 CellChangeTimes->{
  3.8227124476275525`*^9, {3.8227127107271605`*^9, 3.822712744343157*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData["f"], "Input",
 CellChangeTimes->{3.8227127830543947`*^9}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{
    RowBox[{"-", "n"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "+", 
      SuperscriptBox["n", "2"]}], ")"}]}], "+", 
   RowBox[{"2", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "-", 
      SuperscriptBox["n", "2"]}], ")"}], " ", "p", " ", 
    RowBox[{"ArcTan", "[", 
     FractionBox[
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{"k", "+", "p"}], ")"}]}], "n"], "]"}]}], "-", 
   RowBox[{"4", " ", "k", " ", "n", " ", "p", " ", 
    RowBox[{"Log", "[", 
     FractionBox[
      SuperscriptBox["p", "2"], 
      RowBox[{
       SuperscriptBox["n", "2"], "+", 
       RowBox[{"4", " ", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{"k", "+", "p"}], ")"}], "2"]}]}]], "]"}]}]}], 
  RowBox[{
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "+", 
      SuperscriptBox["n", "2"]}], ")"}], "2"], " ", "p"}]]], "Output",
 CellChangeTimes->{3.822712783391674*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"f", "=", 
  RowBox[{"f", "*", 
   RowBox[{"2", "/", "\[Pi]"}]}]}]], "Input",
 CellChangeTimes->{{3.8227130820854073`*^9, 3.8227130867117395`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{"2", " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "n"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"4", " ", 
         SuperscriptBox["k", "2"]}], "+", 
        SuperscriptBox["n", "2"]}], ")"}]}], "+", 
     RowBox[{"2", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"4", " ", 
         SuperscriptBox["k", "2"]}], "-", 
        SuperscriptBox["n", "2"]}], ")"}], " ", "p", " ", 
      RowBox[{"ArcTan", "[", 
       FractionBox[
        RowBox[{"2", " ", 
         RowBox[{"(", 
          RowBox[{"k", "+", "p"}], ")"}]}], "n"], "]"}]}], "-", 
     RowBox[{"4", " ", "k", " ", "n", " ", "p", " ", 
      RowBox[{"Log", "[", 
       FractionBox[
        SuperscriptBox["p", "2"], 
        RowBox[{
         SuperscriptBox["n", "2"], "+", 
         RowBox[{"4", " ", 
          SuperscriptBox[
           RowBox[{"(", 
            RowBox[{"k", "+", "p"}], ")"}], "2"]}]}]], "]"}]}]}], ")"}]}], 
  RowBox[{
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "+", 
      SuperscriptBox["n", "2"]}], ")"}], "2"], " ", "p", " ", 
   "\[Pi]"}]]], "Output",
 CellChangeTimes->{3.822713087249523*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"D", "[", 
    RowBox[{"f", ",", "p"}], "]"}], "-", 
   RowBox[{
    RowBox[{"an", "[", 
     RowBox[{
      RowBox[{"k", "+", "p"}], ",", "n"}], "]"}], "/", 
    RowBox[{"p", "^", "2"}]}]}], "//", "FullSimplify"}]], "Input",
 CellChangeTimes->{{3.8227127842906046`*^9, 3.822712789514861*^9}, 
   3.8227128215823336`*^9, 3.8227130942223625`*^9}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{3.8227127897891273`*^9, 3.8227128226654067`*^9, 
  3.8227130966367235`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"f", "//", "Simplify"}]], "Input",
 CellChangeTimes->{{3.822713099266853*^9, 3.822713104955678*^9}}],

Cell[BoxData[
 RowBox[{"-", 
  FractionBox[
   RowBox[{"2", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"-", "4"}], " ", 
          SuperscriptBox["k", "2"]}], "+", 
         SuperscriptBox["n", "2"]}], ")"}], " ", "p", " ", 
       RowBox[{"ArcTan", "[", 
        FractionBox[
         RowBox[{"2", " ", 
          RowBox[{"(", 
           RowBox[{"k", "+", "p"}], ")"}]}], "n"], "]"}]}], "+", 
      RowBox[{"n", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"4", " ", 
          SuperscriptBox["k", "2"]}], "+", 
         SuperscriptBox["n", "2"], "+", 
         RowBox[{"4", " ", "k", " ", "p", " ", 
          RowBox[{"Log", "[", 
           FractionBox[
            SuperscriptBox["p", "2"], 
            RowBox[{
             SuperscriptBox["n", "2"], "+", 
             RowBox[{"4", " ", 
              SuperscriptBox[
               RowBox[{"(", 
                RowBox[{"k", "+", "p"}], ")"}], "2"]}]}]], "]"}]}]}], 
        ")"}]}]}], ")"}]}], 
   RowBox[{
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       RowBox[{"4", " ", 
        SuperscriptBox["k", "2"]}], "+", 
       SuperscriptBox["n", "2"]}], ")"}], "2"], " ", "p", " ", 
    "\[Pi]"}]]}]], "Output",
 CellChangeTimes->{{3.822713102397172*^9, 3.82271310568087*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"fint", "[", 
    RowBox[{"k", ",", "p", ",", "n"}], "]"}], "+", "f"}], "//", 
  "Simplify"}]], "Input",
 CellChangeTimes->{{3.822713153690051*^9, 3.822713171323168*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{"8", " ", 
    SuperscriptBox["k", "2"]}], "-", 
   RowBox[{"2", " ", 
    SuperscriptBox["n", "2"]}]}], 
  SuperscriptBox[
   RowBox[{"(", 
    RowBox[{
     RowBox[{"4", " ", 
      SuperscriptBox["k", "2"]}], "+", 
     SuperscriptBox["n", "2"]}], ")"}], "2"]]], "Output",
 CellChangeTimes->{{3.822713162617228*^9, 3.8227131717818894`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"an", "[", 
  RowBox[{"x", ",", "n"}], "]"}]], "Input",
 CellChangeTimes->{{3.822712830296524*^9, 3.822712832110276*^9}}],

Cell[BoxData[
 FractionBox["n", 
  RowBox[{"2", " ", "\[Pi]", " ", 
   RowBox[{"(", 
    RowBox[{
     FractionBox[
      SuperscriptBox["n", "2"], "4"], "+", 
     SuperscriptBox["x", "2"]}], ")"}]}]]], "Output",
 CellChangeTimes->{3.822712832423438*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData["f"], "Input",
 CellChangeTimes->{3.822712842375678*^9}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{
    RowBox[{"-", "n"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "+", 
      SuperscriptBox["n", "2"]}], ")"}]}], "+", 
   RowBox[{"2", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "-", 
      SuperscriptBox["n", "2"]}], ")"}], " ", "p", " ", 
    RowBox[{"ArcTan", "[", 
     FractionBox[
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{"k", "+", "p"}], ")"}]}], "n"], "]"}]}], "-", 
   RowBox[{"4", " ", "k", " ", "n", " ", "p", " ", 
    RowBox[{"Log", "[", 
     FractionBox[
      SuperscriptBox["p", "2"], 
      RowBox[{
       SuperscriptBox["n", "2"], "+", 
       RowBox[{"4", " ", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{"k", "+", "p"}], ")"}], "2"]}]}]], "]"}]}]}], 
  RowBox[{
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "+", 
      SuperscriptBox["n", "2"]}], ")"}], "2"], " ", "p"}]]], "Output",
 CellChangeTimes->{3.8227128427845864`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"f", "/.", 
  RowBox[{"p", "\[Rule]", "Infinity"}]}]], "Input",
 CellChangeTimes->{{3.8227128443144045`*^9, 3.822712850525908*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Infinity", "::", "indet"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Indeterminate expression \[NoBreak]\\!\\(0\\\\ \[Infinity]\
\\)\[NoBreak] encountered. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/indet\\\", ButtonNote -> \
\\\"Infinity::indet\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.822712852941246*^9}],

Cell[BoxData["Indeterminate"], "Output",
 CellChangeTimes->{3.822712852942243*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData["f"], "Input",
 CellChangeTimes->{3.8227128618926754`*^9}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{
    RowBox[{"-", "n"}], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "+", 
      SuperscriptBox["n", "2"]}], ")"}]}], "+", 
   RowBox[{"2", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "-", 
      SuperscriptBox["n", "2"]}], ")"}], " ", "p", " ", 
    RowBox[{"ArcTan", "[", 
     FractionBox[
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{"k", "+", "p"}], ")"}]}], "n"], "]"}]}], "-", 
   RowBox[{"4", " ", "k", " ", "n", " ", "p", " ", 
    RowBox[{"Log", "[", 
     FractionBox[
      SuperscriptBox["p", "2"], 
      RowBox[{
       SuperscriptBox["n", "2"], "+", 
       RowBox[{"4", " ", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{"k", "+", "p"}], ")"}], "2"]}]}]], "]"}]}]}], 
  RowBox[{
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "+", 
      SuperscriptBox["n", "2"]}], ")"}], "2"], " ", "p"}]]], "Output",
 CellChangeTimes->{3.8227128622098694`*^9}]
}, Open  ]],

Cell[BoxData["f"], "Input",
 CellChangeTimes->{3.8227132750003796`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"ff", "[", 
   RowBox[{"k_", ",", "p_", ",", "n_"}], "]"}], ":=", 
  FractionBox[
   RowBox[{"2", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"-", "n"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"4", " ", 
          SuperscriptBox["k", "2"]}], "+", 
         SuperscriptBox["n", "2"]}], ")"}]}], "+", 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"4", " ", 
          SuperscriptBox["k", "2"]}], "-", 
         SuperscriptBox["n", "2"]}], ")"}], " ", "p", " ", 
       RowBox[{"ArcTan", "[", 
        FractionBox[
         RowBox[{"2", " ", 
          RowBox[{"(", 
           RowBox[{"k", "+", "p"}], ")"}]}], "n"], "]"}]}], "-", 
      RowBox[{"4", " ", "k", " ", "n", " ", "p", " ", 
       RowBox[{"Log", "[", 
        FractionBox[
         SuperscriptBox["p", "2"], 
         RowBox[{
          SuperscriptBox["n", "2"], "+", 
          RowBox[{"4", " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"k", "+", "p"}], ")"}], "2"]}]}]], "]"}]}]}], ")"}]}], 
   RowBox[{
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       RowBox[{"4", " ", 
        SuperscriptBox["k", "2"]}], "+", 
       SuperscriptBox["n", "2"]}], ")"}], "2"], " ", "p", " ", 
    "\[Pi]"}]]}]], "Input",
 CellChangeTimes->{{3.8227133989629173`*^9, 3.8227134148780174`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"D", "[", 
    RowBox[{"f", ",", "p"}], "]"}], "-", 
   RowBox[{
    RowBox[{"an", "[", 
     RowBox[{
      RowBox[{"k", "+", "p"}], ",", "n"}], "]"}], "/", 
    RowBox[{"p", "^", "2"}]}]}], "//", "FullSimplify"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ff", "[", 
  RowBox[{"1.0", ",", "10000000", ",", "4"}], "]"}]], "Input",
 CellChangeTimes->{{3.8227134190728087`*^9, 3.8227134927825165`*^9}, {
  3.822713529521511*^9, 3.8227135296870627`*^9}}],

Cell[BoxData[
 RowBox[{"-", "0.02469830397557575`"}]], "Output",
 CellChangeTimes->{
  3.822713287685214*^9, {3.8227134525509443`*^9, 3.8227134933269997`*^9}, 
   3.8227135305990376`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"an", "[", 
  RowBox[{
   RowBox[{"1", "+", "1"}], ",", "4"}], "]"}]], "Input",
 CellChangeTimes->{{3.8227134533803463`*^9, 3.822713458821164*^9}}],

Cell[BoxData[
 FractionBox["1", 
  RowBox[{"4", " ", "\[Pi]"}]]], "Output",
 CellChangeTimes->{3.822713459858328*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"N", "[", "%", "]"}]], "Input",
 CellChangeTimes->{{3.8227134607599163`*^9, 3.822713461927455*^9}}],

Cell[BoxData["0.07957747154594767`"], "Output",
 CellChangeTimes->{3.8227134621827717`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"an", "[", 
   RowBox[{
    RowBox[{"k", "+", "p"}], ",", "n"}], "]"}], "/", 
  RowBox[{"p", "^", "2"}]}]], "Input",
 CellChangeTimes->{{3.8227132902086835`*^9, 3.8227132958651276`*^9}}],

Cell[BoxData[
 FractionBox["n", 
  RowBox[{"2", " ", 
   SuperscriptBox["p", "2"], " ", 
   RowBox[{"(", 
    RowBox[{
     FractionBox[
      SuperscriptBox["n", "2"], "4"], "+", 
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{"k", "+", "p"}], ")"}], "2"]}], ")"}], " ", 
   "\[Pi]"}]]], "Output",
 CellChangeTimes->{3.8227132962511554`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"fint", "[", 
   RowBox[{"k_", ",", "kc_", ",", "n_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"n", "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"4", "*", 
            RowBox[{"k", "^", "2"}]}], "+", 
           RowBox[{"n", "^", "2"}]}], ")"}]}], "-", 
        RowBox[{"2", "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"4", "*", 
            RowBox[{"k", "^", "2"}]}], "-", 
           RowBox[{"n", "^", "2"}]}], ")"}], "*", "kc", "*", 
         RowBox[{"ArcTan", "[", 
          RowBox[{
           RowBox[{"2", "/", "n"}], "*", 
           RowBox[{"(", 
            RowBox[{"k", "+", "kc"}], ")"}]}], "]"}]}], "+", 
        RowBox[{"4", "*", "n", "*", "k", "*", "kc", "*", 
         RowBox[{"Log", "[", 
          RowBox[{
           RowBox[{"kc", "^", "2"}], "/", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"n", "^", "2"}], "+", 
             RowBox[{"4", "*", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{"k", "+", "kc"}], ")"}], "^", "2"}]}]}], ")"}]}], 
          "]"}]}]}], ")"}], "/", "kc"}], "/", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"n", "^", "2"}], "+", 
        RowBox[{"4", "*", 
         RowBox[{"k", "^", "2"}]}]}], ")"}], "^", "2"}]}], "*", 
    RowBox[{"2", "/", "\[Pi]"}]}], "+", 
   RowBox[{"2", "*", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"4", "*", 
        RowBox[{"k", "^", "2"}]}], "-", 
       RowBox[{"n", "^", "2"}]}], ")"}], "/", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"4", "*", 
         RowBox[{"k", "^", "2"}]}], "+", 
        RowBox[{"n", "^", "2"}]}], ")"}], "^", "2"}]}]}], "+", 
   RowBox[{"4", "*", "k", "*", 
    RowBox[{"n", "/", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"4", "*", 
         RowBox[{"k", "^", "2"}]}], "+", 
        RowBox[{"n", "^", "2"}]}], ")"}], "^", "2"}]}], "*", 
    RowBox[{"Log", "[", "4", "]"}], "*", 
    RowBox[{"2", "/", "\[Pi]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.8227135920872583`*^9, 3.822713608978645*^9}, {
  3.8227137619488635`*^9, 3.8227137636806355`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"fint", "[", 
  RowBox[{"k", ",", 
   SubscriptBox["k", "c"], ",", "n"}], "]"}]], "Input",
 CellChangeTimes->{{3.822713716030431*^9, 3.822713723352302*^9}}],

Cell[BoxData[
 RowBox[{
  FractionBox[
   RowBox[{"2", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "-", 
      SuperscriptBox["n", "2"]}], ")"}]}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "+", 
      SuperscriptBox["n", "2"]}], ")"}], "2"]], "+", 
  FractionBox[
   RowBox[{"8", " ", "k", " ", "n", " ", 
    RowBox[{"Log", "[", "4", "]"}]}], 
   RowBox[{
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       RowBox[{"4", " ", 
        SuperscriptBox["k", "2"]}], "+", 
       SuperscriptBox["n", "2"]}], ")"}], "2"], " ", "\[Pi]"}]], "+", 
  FractionBox[
   RowBox[{"2", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"n", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"4", " ", 
          SuperscriptBox["k", "2"]}], "+", 
         SuperscriptBox["n", "2"]}], ")"}]}], "-", 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"4", " ", 
          SuperscriptBox["k", "2"]}], "-", 
         SuperscriptBox["n", "2"]}], ")"}], " ", 
       RowBox[{"ArcTan", "[", 
        FractionBox[
         RowBox[{"2", " ", 
          RowBox[{"(", 
           RowBox[{"k", "+", 
            SubscriptBox["k", "c"]}], ")"}]}], "n"], "]"}], " ", 
       SubscriptBox["k", "c"]}], "+", 
      RowBox[{"4", " ", "k", " ", "n", " ", 
       RowBox[{"Log", "[", 
        FractionBox[
         SubsuperscriptBox["k", "c", "2"], 
         RowBox[{
          SuperscriptBox["n", "2"], "+", 
          RowBox[{"4", " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"k", "+", 
              SubscriptBox["k", "c"]}], ")"}], "2"]}]}]], "]"}], " ", 
       SubscriptBox["k", "c"]}]}], ")"}]}], 
   RowBox[{
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       RowBox[{"4", " ", 
        SuperscriptBox["k", "2"]}], "+", 
       SuperscriptBox["n", "2"]}], ")"}], "2"], " ", "\[Pi]", " ", 
    SubscriptBox["k", "c"]}]]}]], "Output",
 CellChangeTimes->{3.822713724102071*^9, 3.8227137671026297`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"fint", "[", 
   RowBox[{"1", ",", "2", ",", "4"}], "]"}], "//", "N"}]], "Input",
 CellChangeTimes->{{3.8227136101881266`*^9, 3.822713618095186*^9}, {
  3.8227136894025717`*^9, 3.822713690430465*^9}}],

Cell[BoxData["0.011187729200601043`"], "Output",
 CellChangeTimes->{3.8227136185768967`*^9, 3.8227136913931117`*^9, 
  3.822713769966854*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"NIntegrate", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"an", "[", 
     RowBox[{
      RowBox[{"1", "+", "p"}], ",", "4"}], "]"}], "/", 
    RowBox[{"p", "^", "2"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"p", ",", "2", ",", "Infinity"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8227136280974007`*^9, 3.8227136726295505`*^9}}],

Cell[BoxData["0.011187729200627137`"], "Output",
 CellChangeTimes->{{3.8227136546988974`*^9, 3.822713673335683*^9}, 
   3.8227137764370303`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"k", "=", 
   RowBox[{"-", 
    RowBox[{"RandomReal", "[", "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"kc", "=", 
   RowBox[{"RandomReal", "[", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"n", "=", 
   RowBox[{"RandomReal", "[", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Print", "[", 
   RowBox[{"fint", "[", 
    RowBox[{"k", ",", "kc", ",", "n"}], "]"}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  RowBox[{"NIntegrate", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"an", "[", 
      RowBox[{
       RowBox[{"k", "+", "p"}], ",", "n"}], "]"}], "/", 
     RowBox[{"p", "^", "2"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"p", ",", "kc", ",", "Infinity"}], "}"}]}], "]"}], 
  "]"}]}], "Input",
 CellChangeTimes->{{3.822713788429907*^9, 3.822713845444709*^9}, 
   3.8227139157862015`*^9}],

Cell[CellGroupData[{

Cell[BoxData["0.7518733286065291`"], "Print",
 CellChangeTimes->{{3.822713845846566*^9, 3.8227138601494975`*^9}, {
   3.8227139163536525`*^9, 3.822713919307671*^9}, 3.8227274783617086`*^9}],

Cell[BoxData["0.751873328606566`"], "Print",
 CellChangeTimes->{{3.822713845846566*^9, 3.8227138601494975`*^9}, {
   3.8227139163536525`*^9, 3.822713919307671*^9}, 3.8227274783646927`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Clear", "[", "k", "]"}], ";", 
  RowBox[{"Clear", "[", "n", "]"}], ";", 
  RowBox[{"Clear", "[", "kc", "]"}], ";", 
  RowBox[{"fint", "[", 
   RowBox[{"k", ",", 
    SubscriptBox["k", "c"], ",", "n"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8227274816516438`*^9, 3.82272754512259*^9}, {
  3.822727670144861*^9, 3.8227276750469913`*^9}}],

Cell[BoxData[
 RowBox[{
  FractionBox[
   RowBox[{"2", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "-", 
      SuperscriptBox["n", "2"]}], ")"}]}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"4", " ", 
       SuperscriptBox["k", "2"]}], "+", 
      SuperscriptBox["n", "2"]}], ")"}], "2"]], "+", 
  FractionBox[
   RowBox[{"8", " ", "k", " ", "n", " ", 
    RowBox[{"Log", "[", "4", "]"}]}], 
   RowBox[{
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       RowBox[{"4", " ", 
        SuperscriptBox["k", "2"]}], "+", 
       SuperscriptBox["n", "2"]}], ")"}], "2"], " ", "\[Pi]"}]], "+", 
  FractionBox[
   RowBox[{"2", " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"n", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"4", " ", 
          SuperscriptBox["k", "2"]}], "+", 
         SuperscriptBox["n", "2"]}], ")"}]}], "-", 
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"4", " ", 
          SuperscriptBox["k", "2"]}], "-", 
         SuperscriptBox["n", "2"]}], ")"}], " ", 
       RowBox[{"ArcTan", "[", 
        FractionBox[
         RowBox[{"2", " ", 
          RowBox[{"(", 
           RowBox[{"k", "+", 
            SubscriptBox["k", "c"]}], ")"}]}], "n"], "]"}], " ", 
       SubscriptBox["k", "c"]}], "+", 
      RowBox[{"4", " ", "k", " ", "n", " ", 
       RowBox[{"Log", "[", 
        FractionBox[
         SubsuperscriptBox["k", "c", "2"], 
         RowBox[{
          SuperscriptBox["n", "2"], "+", 
          RowBox[{"4", " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"k", "+", 
              SubscriptBox["k", "c"]}], ")"}], "2"]}]}]], "]"}], " ", 
       SubscriptBox["k", "c"]}]}], ")"}]}], 
   RowBox[{
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       RowBox[{"4", " ", 
        SuperscriptBox["k", "2"]}], "+", 
       SuperscriptBox["n", "2"]}], ")"}], "2"], " ", "\[Pi]", " ", 
    SubscriptBox["k", "c"]}]]}]], "Output",
 CellChangeTimes->{
  3.822727491363246*^9, {3.8227275340536275`*^9, 3.822727545929369*^9}, 
   3.822727675771758*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"fint", "[", 
    RowBox[{"2", ",", "kc", ",", "4"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"kc", ",", "1", ",", "20"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8227275769799824`*^9, 3.8227275933481784`*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
    1.], LineBox[CompressedData["
1:eJwV1Hs41NkfB/CvyRpyn3EbY0VK0tJILkXz+UTjMooJyWUSswmpnWnaSqUa
JbRdtnRTblO2CSXVYLMkhe0i0kb8WokuSlEq2qHod/rj+3yf1/P+nM/znHOe
87EWiYPjaBRFFZHv+/+80UP7vJ10dG3b31jSlwTUrXyBSyod7dh/5SY+SwJF
F/pwTtIxW5lSPnEmCT7S9zjZ3qHjENvcI8k2CTKjdTUZMzVx6Nv92nGXNVCp
bVnR/1oT6ywVy5QbE2CvyfvZr4Y1sbVld6FWVAIIretKXlBaaPC0ZCIdE4Dm
JpI/NdPCzoasmhHtBFgiUuxr99PCmGkm2Tl/xMOLq46iG0VaKNF43X6vazUY
rgb9kwmTsS41/uik+DjICn+w58RG4lKfju7AOGAE/Pz12K7JaMcrOtvtSszJ
6M/KnYxWa/Na+HTisZb6/a2TEcP7OpqLV4HRwehkmbs2YvT450UjP4NZ5fbe
eE0d7KlpyvT5QwTZRfrhq4110C+1SOJ2SARmOadbVk3VwVbaCh9eighYssbq
WE8dFAzk91eGisCcr3s8SqKDHGXD2gENEbC7c/lBnTooud7016+SWLDUqFG6
Fuli3ZsjyVfCYkDDeWCnqlwXs5PZ1U4+MfBupcXiv27oYqZvc9Atlxi4XpXy
YsG/xIdZnXONYyBm3QJjnp4eyvYEye1LVsKZttpNIb8SB9mtzr4cDbaFN+ZJ
UB/lpsNveN1C0Gv98MOcJcRvjV4Kbgrh81frfz5FEAdtZckUQrgVtmvN5g36
KNt/kMcXCyFhsnfODoU+HiqLD8mgCaHdXbTDR8MAOeZdKunsKDCVP38autYA
ZVZ0Zc7pCPjFzovuJzXAurCGOSv2RUDDJbmjR7IBXrom3u2+MQKkN1Zst04j
9WdC1Dz9I6DlWYf5uxzSL2oiVvoxHNKm31uW0UTqPy2x3O4fDh9Kypuq7A0x
ZoG/p7nWcvBzZn66wDHEIQ/zsgcjYZBfvd5c7mqIh04mVuQ/C4PF9xwT072I
1TImbakJg6LBIo2QSEOUU6Oxs9aHQTQnb+HAXkOUReZeWNe9DJoq9vxp2U/q
q9SWK26FgupygzDrvSHWGSxzEleEgu3FSTSNzyR/FxGzpDAUZIrUxYM0Bh7y
jwH/naHgemJ7b7UFAyXb6057uIfCmeSN2hECBtaFXCsZKg2BrR5xK4/8SfKL
sn8fnguGIrc/1Om1DBQEaB/Syw6GdufnxVsbGDg0z5QRlRkMnJ9ih2MfMNDA
1tx8ZmIw9Fms2Mt5y0CO9tCu47OCIXgiVNlsyUSOe+Te08qlYH9jEV0znYmt
f24w7rgvgPmz7w7q7meiwH+wSqNeAPy8oDZmFnHX3ATvSgEkbYk8PSWfiUjd
KXieK4ALThIPt0omSgKby66tEYDDmRzx6j4mWr17xL+rKQDO7o8dDT5GSPVv
0i/PCQRXnrwolW5M9ufGlH0IgK2tw/MSdY0R/c+G9zwJgNoo/yYB0xitxg+o
hdwNAF/px8EpU4yx57eCGeLCAAgvWORc60rWbzOUyMLI+tHX18bijJEKHCgd
vs6Ha6VODzc0kHzOsSMnCvxhkUn9eFyqCeL8BZ7VB3xhpbFKpyXDBKnXIfW3
U3xhq5GDhetB4tmFB/qSfOES48R8jRwTlDn35wXzfYGtn7RZUW6CVnttDr+h
+5J5xfzQ98oE5d2PPLXSfUA+KnoRH2iKVEXyi10HeDDeRWtKZJthjHKT7Y1G
bxhbsoO9ypo4OV/qUeUNqtrRpOgZxKtuL6274A3D8iGdEGczpMJ25fYf8YY3
q7oDPQPMsKd4oWOVyBseDVT9o7fNDOWsD41eNG+49FXyWPnYjMwXSU6lnxeI
LJ72j59koczhp3jqA4KrmsfVZDkL5dlh5iMvESb3HU//pGBhz98hqz89RlCW
Bdr0K1kYM/vHAp1GBHXvWmFbMwvRa+GGSycRFIl590vUzLHnAf1e+SKEN5VR
FcsTzVGWf2tqpxfAhqWdssvubKQcSq1cXnlCSsikMV1gY10EzaCk3RPSljlu
WMMjeWMwx6HBE45GpMXZBLMRLQ/eXXHaEypjOfxjSWyUTeroWC30BNX6TOaW
fJK7tlnptHtASpa7AtUtsO7sjem81vmw52F20/37Fojdbn6X37qDxHdyem27
BVKXhRe9u90hsnoblv5L6sdVV5+2ugOnMLrit1fEm2gn51a6Q5d0WgHvmwXK
FP5HBanuMJdxSVrt+CPWie2N+Sx3eCG4xTp3gPhh2biWwA28W0bidwRYkvvf
fE2vwwXUK33lg0stETOfaXXedYG/8052CsMtsY4/S1JS6wL8dQv8PeKIqzYn
rFG4gEBnj71qB/GKCYuwTS4QzTcaXH+F+PBCYbupCyT/PUcax5pCzk8CHrFz
4UKtOGXxK+I03iVKxxl6md7r25VWiC0LuYcfc8CrXhhjlGaN1LfehLsXHWAa
K7I4UzAVqXnx65pD7UF4fjqbp2+DVG/2Oq2zM2DN1W9+ht02KBuYmGBYT4cW
taKC6QXTyLw8zasxtYErH5wFYuF0lHUGh5ctt4a1wpaa2zNtkdIvtumKnAIp
k9yC1/QT65Xjl3wLEE2/WMiumkGee9b+ZWrmcGdQM4i2xw6pp0/OfKkxhfSP
TjZO3jOR2vCkb+CxMbCkiUqxoT1Sa//XbahpBDRH8drPD4jF4vrxVAYkHQrf
eeXgLKQcj+gYzDSEknWCwB1RPyG138lzQZk+QIaJT5a5A1KnSm22bdKFEvbb
0Yp+4ppDI3F+2pA0cOCR9XlHpIRrpX1iLSj92Ui3PmE2Una9FtPu0OFZSdu3
eWwOUpa/tL231QCWV/oB123ElfVhW26qQ7UJjM5/TvzIiWMZOwnyavsNI/2d
kMp8kq4xjQYlaV+4eUXEp7IeWxupwdRMI/5r/TlIKX78vVGdgsbio4nLJcSa
zcXnXCe4Rfkarr0txDV9nyoPfuV+bt3ZI3VxRmpR7z7myzGu2DqmwO4YsSxp
I2fjKDc05/bXh8eJ7Q6PqYlHufONnSJ3ZhN73tzZnjDKVdeiGT/KIaa0M3YJ
R7mnhgp/232GePi8dNh7lNtw/dXG7jLiHlNtN+Yo1yxavPjoXeJfLWqclCru
eEdHMdwjvs2uWliq4j5bivS3zd/7/ZAXck7FLeUZ3lz4gLhmr8n2HBXXy0Hp
9r6DOM2mum+3irvu6+ep/JfEwnFro2UqbvCmlbKRPmJNXWevIBXXfejWE/lr
Yh31KVJ/FZf2Ijv7v7ffXZ7WsUDFzW6ar3f2I/HVXdKrtirudl5hkmCYONfJ
Z8RKxRVd177zZYS47fU7Z7aK66Ds2h08+r1/Rk2lvorLcOD1jo99P6+xp2Na
Kq5KUcot/kqs/ns7qqu4T6xMckMniC/sy82c+I9789SO0W/fiIfc5v6j+o/7
fy/ybyc=
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{5., 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  Method->{
   "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None},
  PlotRange->{{1, 20}, {0., 0.0034058890709283375`}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.822727594282499*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Refine", "[", 
  RowBox[{
   RowBox[{"Series", "[", 
    RowBox[{
     RowBox[{"fint", "[", 
      RowBox[{"k", ",", "kc", ",", "n"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"kc", ",", "Infinity", ",", "8"}], "}"}]}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"k", "\[Element]", "Reals"}], ",", 
     RowBox[{"n", ">", "0"}]}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8227276385610266`*^9, 3.822727662015358*^9}, {
  3.8227277497544203`*^9, 3.8227277856468687`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
   FractionBox["n", 
    RowBox[{"6", " ", "\[Pi]", " ", 
     SuperscriptBox["kc", "3"]}]], "-", 
   FractionBox[
    RowBox[{"k", " ", "n"}], 
    RowBox[{"4", " ", "\[Pi]", " ", 
     SuperscriptBox["kc", "4"]}]], "+", 
   FractionBox[
    RowBox[{
     RowBox[{"12", " ", 
      SuperscriptBox["k", "2"], " ", "n"}], "-", 
     SuperscriptBox["n", "3"]}], 
    RowBox[{"40", " ", "\[Pi]", " ", 
     SuperscriptBox["kc", "5"]}]], "+", 
   FractionBox[
    RowBox[{
     RowBox[{
      RowBox[{"-", "4"}], " ", 
      SuperscriptBox["k", "3"], " ", "n"}], "+", 
     RowBox[{"k", " ", 
      SuperscriptBox["n", "3"]}]}], 
    RowBox[{"12", " ", "\[Pi]", " ", 
     SuperscriptBox["kc", "6"]}]], "+", 
   FractionBox[
    RowBox[{
     RowBox[{"80", " ", 
      SuperscriptBox["k", "4"], " ", "n"}], "-", 
     RowBox[{"40", " ", 
      SuperscriptBox["k", "2"], " ", 
      SuperscriptBox["n", "3"]}], "+", 
     SuperscriptBox["n", "5"]}], 
    RowBox[{"224", " ", "\[Pi]", " ", 
     SuperscriptBox["kc", "7"]}]], "+", 
   FractionBox[
    RowBox[{
     RowBox[{
      RowBox[{"-", "48"}], " ", 
      SuperscriptBox["k", "5"], " ", "n"}], "+", 
     RowBox[{"40", " ", 
      SuperscriptBox["k", "3"], " ", 
      SuperscriptBox["n", "3"]}], "-", 
     RowBox[{"3", " ", "k", " ", 
      SuperscriptBox["n", "5"]}]}], 
    RowBox[{"128", " ", "\[Pi]", " ", 
     SuperscriptBox["kc", "8"]}]], "+", 
   InterpretationBox[
    SuperscriptBox[
     RowBox[{"O", "[", 
      FractionBox["1", "kc"], "]"}], "9"],
    SeriesData[$CellContext`kc, 
     DirectedInfinity[1], {}, 3, 9, 1],
    Editable->False]}],
  SeriesData[$CellContext`kc, 
   DirectedInfinity[1], {
   Rational[1, 6] $CellContext`n/Pi, 
    Rational[-1, 4] $CellContext`k $CellContext`n/Pi, 
    Rational[1, 40] (12 $CellContext`k^2 $CellContext`n - $CellContext`n^3)/
    Pi, Rational[
     1, 12] ((-4) $CellContext`k^3 $CellContext`n + $CellContext`k \
$CellContext`n^3)/Pi, 
    Rational[1, 224] (80 $CellContext`k^4 $CellContext`n - 
     40 $CellContext`k^2 $CellContext`n^3 + $CellContext`n^5)/Pi, 
    Rational[1, 128] ((-48) $CellContext`k^5 $CellContext`n + 
     40 $CellContext`k^3 $CellContext`n^3 - 
     3 $CellContext`k $CellContext`n^5)/Pi}, 3, 9, 1],
  Editable->False]], "Output",
 CellChangeTimes->{{3.82272766276254*^9, 3.8227276870750694`*^9}, {
  3.822727775892643*^9, 3.822727792499372*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Series", "[", 
  RowBox[{
   RowBox[{"an", "[", 
    RowBox[{"x", ",", "n"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "Infinity", ",", "4"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8227278344862385`*^9, 3.8227278564401517`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
   FractionBox["n", 
    RowBox[{"2", " ", "\[Pi]", " ", 
     SuperscriptBox["x", "2"]}]], "-", 
   FractionBox[
    SuperscriptBox["n", "3"], 
    RowBox[{"8", " ", "\[Pi]", " ", 
     SuperscriptBox["x", "4"]}]], "+", 
   InterpretationBox[
    SuperscriptBox[
     RowBox[{"O", "[", 
      FractionBox["1", "x"], "]"}], "5"],
    SeriesData[$CellContext`x, 
     DirectedInfinity[1], {}, 2, 5, 1],
    Editable->False]}],
  SeriesData[$CellContext`x, 
   DirectedInfinity[1], {
   Rational[1, 2] $CellContext`n/Pi, 0, Rational[-1, 8] $CellContext`n^3/Pi}, 
   2, 5, 1],
  Editable->False]], "Output",
 CellChangeTimes->{3.822727857084945*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Refine", "[", 
  RowBox[{
   RowBox[{"Integrate", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"an", "[", 
       RowBox[{
        RowBox[{"k", "-", "p"}], ",", "1"}], "]"}], "*", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"p", "^", "2"}], "-", "1"}], ")"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"p", ",", 
       RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"k", "\[Element]", "Reals"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.822727978284721*^9, 3.8227280266377993`*^9}, {
  3.8227280643867025`*^9, 3.822728085106137*^9}, {3.8227281752757883`*^9, 
  3.8227281817018456`*^9}}],

Cell[BoxData["$Aborted"], "Output",
 CellChangeTimes->{3.822728012809463*^9, 3.822728051537185*^9, 
  3.822728136924449*^9, 3.8227282060038724`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FourierTransform", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{"1", "-", 
      RowBox[{"Sign", "[", 
       RowBox[{
        RowBox[{"x", "^", "2"}], "-", "1"}], "]"}]}], ")"}], "*", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"x", "^", "2"}], "-", "1"}], ")"}]}], ",", "x", ",", "k"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.8227282102432404`*^9, 3.822728259029563*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{"8", " ", "k", " ", 
    RowBox[{"Cos", "[", "k", "]"}]}], "-", 
   RowBox[{"8", " ", 
    RowBox[{"Sin", "[", "k", "]"}]}]}], 
  RowBox[{
   SuperscriptBox["k", "3"], " ", 
   SqrtBox[
    RowBox[{"2", " ", "\[Pi]"}]]}]]], "Output",
 CellChangeTimes->{3.8227282621568184`*^9}]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.822728375763379*^9, 3.822728439320613*^9}, 
   3.8227284822793612`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{"1", "-", 
      RowBox[{"Sign", "[", 
       RowBox[{
        RowBox[{"x", "^", "2"}], "-", "1"}], "]"}]}], ")"}], "*", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"x", "^", "2"}], "-", "1"}], ")"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", 
     RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.822728328670276*^9, 3.8227283360354986`*^9}}],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
     1.], LineBox[CompressedData["
1:eJxFx2kogwEcx/HHQiYzIUdNrnKVITVXLEfu+7Zn2cgRKwpvSIgUyU2mZtIc
KwvlSJIRJdcr5Sq3tBk1NufKkfj/f/Xt18cxvzy1kEIQBPOn32db3X190VeJ
v6WXura8f6Br+v1dNM/oTRUn/1aF5guHT3YP0d2PbjuD02itJEjmy0MvGeWV
Fa+Yguv3DlRreTSwcszW2o9hAvY4j6bST43BoxYS7XEjFbxObrfL2UZgxjDD
QXhpCI4LJ6v3Ww3A9rqiBvcYfbDn1FBDxQMFfFE1wjyQ6oFLRN0SUSkBzq7L
cZY5fcr/7aZyJ9PFOrAgMbn56voNzKlMCmaxXsC9ckWyfEEDFtcSNOUkWhpg
s2U+gl6ZjQotakMrJsZ9aHx0SFeBOUlFK/kX+y/kE1hr99bpkIL+PKHHx0ai
LTPYG0Pe6JAY8VyEAbrPh9vfM/UIZlOOMpl66FvXXIXZqxrckXBTrblHnwnV
osUjdJ0X9TpsBr3MDSzP4qELGleJwAy0iTSyhxGH5mpS5y9ZaF2rQCcwRTcN
pKkrjNHfZ2vINQ==
      "]], LineBox[CompressedData["
1:eJw1l3k01V3Ux1EkZahHkgZDpCc9pZSh6KgkRAMqKSpF9CRDIaIIRSk9Kg2G
MhaRqYzxNVWGKyVzXDJkuu5Fpntxf+9913rfP84667P23mud/d17n7OOvI2T
qa0AHx9fL2/97z7zzUJWZogF38/LS0Q7AvB7Yua0cicLfIG0A18qQrF9t7A/
o5rHd6zjHow/R0CIZGJGDo/9TQ5voGJQ3yRX6R7Hi++o+drK9xpXL2uJ819n
QfVavpKpTzrK8/Zt/WzHwpm6mnhrdiYkBU2P3jPl2XfWdX+/8gFZERcjlm1g
oVPzXQzjXD7GKiKUN7Qw0exy4oGndCmc1nKJmQYTfN5TCxRQjYvNpc4JfxiI
v03z8tFqguOPs2V57QxItomdGWxqgvNXPqmvXxi49MldscOtGR7luwomIxj4
ZpZjdCa9BQHp+YIGegwEnZ5ZXCLbhuigzBeDj4dw8Pn+bdwrnfihGVu+WX0Q
yd1jOv+I96JRbfdyPblBMI2HY38d60Xrpk4HC5FBKPh6TZ2I7kWn4hpxP/oA
nPVMnMI3/gZT4sWJutsD8NKrf75Lvw/CA/8xrzb1w17BvvXA5QFoP/eTzr/W
hwdHXy+0/G8YJpahT7ed60OiU1n62W/DsF4ZJZVm0gfrZH3XzeJM3IzKlYxX
6EOm3y291SFMlMawJB7QfoOeU6u0x5+F/clWC23kfmODfKejRdYITAu0OAsr
elDb5y2yxG8M57z3ewZk9qCb6x0wGTmGqzpHp7mRPTBITF+E3DE8gfPkH5ce
CCRm7mazxtBc9nqsfWUPdjjIGy+1+gMr2rKhDKdu5HKTWhu3jMO+ffTnieVd
6LG9QE/6NIEspzZB/XldOJWW2n24bQJc/i+bt7J+4aWKXg59bAKP10X6i3z5
hYCPZyvfyU6i2ElfpcD9F6rI6edy1yYhLRDhtaqxE/tp5/qNFKfw0CRoa9Hf
HZg/F7Yn32Ea61ytbofzd6DiZceSnT7T+Bi+tfVyCx3H3LzTYh9Oo7+jzU82
mI7vuYlH1+VMQ9d1y3ff/nboxRSaUAJsjIS3Ou193QaZ6XHG1BM2DnVuTK1Q
bMW9wuC+eUkcGJ9XPkPYLYhRQG9hPgdG/fJ/Zde0QNvroL4ljQO9ESnPOPcW
CPBf37SDxYEmH/++GxXNqL60KCJcbQYKcg1tapeacLxwYW9J9gxk42tDk3Wb
UGxEF479PINV66v2yC9rwruqxaHnG2ewXBVvxIoasSLWo/ThxAwW6yZdHRBv
xKa0gE+BW2cxcdpn8cusepju8dbuSpzFWLcHlgXVQ2yaFdj/fhYjF1xdQ07V
49Crwi3fS2cx6HSh2VOwHmcr7NIN6bOg3zwSb378B0rdRwyb/prj6aakLTLz
HTeMfw2f95yDOiKaBgO+ovWnJttwExcyi+cPRKt9RUjF8Y1MDS7mLC5xTLtq
8FjN7MG13VyUj2mvLiA1iAvooBNzLszW0W3ucapRuvRg5rNrXDjfl2OqOFdi
fscthSsfuTBvDaI61lRicYbDXq9yLjSURyUe11TgIv8/jy7QuOAWl6jNbqiA
8mpDjP/kImTcxpPW+xmvdJsMW9hcJJ1MmO94shxPNA6sr1alEPJGVEpepBzj
TLb6dg0KzhNuyg25ZXijGHg1UIeCZqi+kY5UGZYs2uLdbEjhc2lfqOj3EvSU
0QsNzlJIFj8cU3KzBExHl709dhTun8rNdNtUApPwzinbSxTMJ4Ma2u8VQ9TG
fHilBy/ezFz91P0ihPT+J2J6j+c/rHiifLAQe24JL5B9yPO/M3F9o0EhjKNO
bq17TKE7P7xkVuAjbA4+XTYURYFPocU46loeHkSeF+t7R+FLQZKTYGMuWCrw
682keHPvFeaolosErcm1VdkUVgWvbNZhZmPJN/b4viIKWiNW5+jn3mOw/1jl
AhqFS/nxAx7lWYi54jSsXEshOmDIaYlSFoRGLu7eXEdh3gpPb73fGWgqFPH8
00RBvRsCdP0M6ImnnUpvpWCfKhTs8Tod1qbqzubtFJ77r1c5+v0dIlKH15l2
UeAuz9YICktF4afVVW97KNik7NUrMEvB/VrVKMZvChsbrK3kG5OwPZeqWztE
4T8HhoP50zcQfVqwXWaYwuScp/sdi9d4dMGyYopJoXjdk9Dh1ngktuv62I9R
UCpQiJSLjMM2/6CXM38o3D2U/sbMKhZOMlmjbhMUWN06H26viYFFXPGVpkkK
r0fC/a/gJZirM5Tkpnn1XpEZ2bI7ClvvBiw4zKawY0/NB1L+Amv7d6y4wKEw
/9/+rwn6z1CiUW9pO0Ph66N5/Ysqn0DU8zDNaJbCs49r+F0PPIJQSprj8jle
fr1aMs01D5H6bXRXNY+Pr0q2834VAv4+KV07LoVBx85nBgl3IMBa49LP4xuQ
qpZMvoWM/oV1RygKCgovOCVW1yFW32QTy2PDT0KnLS+4YHlG0NpWHhfH1IhW
l1mh9uYqSTaPDeKo4kQVPaKm+2gbxeOAWLu+3Cl7sm+MEcDg8cp3+uJ6ru5E
MPxvwRIeewk6Nhzd5ktcVQw/ePP4d5Z6XPDmQBKSeeCJLI/NbPhcCjfcJeYq
WxKSeOcrlqjaNaoUSiqfsHtkeBzeYXg7TCKMMEbjza/w8qNlb3Lc+OAx+UjU
xjN5egg8+Mv886KnROtmfFkLTy9N2+kdZ4Kfk/Pp06V9PD2dtNvlOUKRROfH
5rE2nt6Jf5UKPw6IJqV9hkfyePW47PZjSsrvFZkY1u/0nqIQb5xLT+fGkvO1
A/fzximQ3A11Pjrx5FVyaMpWXr1b10Z9MvJOIEEe0jNhoxQkZvxSetivic3v
OuWFvH55azf+MkMziZgHzReSZVDQr7N7dMMjmYytXCG5apCCT9IBL+mJFKIp
PYTqXgrSy4oce9XeEbrvG6sb3RTe+6qezXRNI2o/DVSkf1GYeCuiuG46nXSf
33qMv42Cij9LpNgti9jR3gwmfaNwaudfwqkK70m29l1h0xre/fJHXfBF7Xsi
U/V0rK2SAvPcDerK39lERltDJaKUQqae6Lhyey7het/ZlMeb3+5Z1VHJu3mk
r6Ey1YM335IfzJn8Gvmk+4AmWzqZgptSZP/PhwVEMvT5+LwYnr5CKu2hekXk
pKujz/X7FByKDrZ6jxaR9iwJ68NBFF64uzY5RIO4LO0UEfbn9duftn5572JS
v25epew1Ch3nFLtnGkrIjRTLppEzFGL1sprSgsuJrcoh7RX/UMiX23dsXWc5
kd71vkNmHYUfs431keqfCN1x3yFqDW8+PnC+B3d/InGKjDBnCQq2Snto53S+
kKUso5TxUS7WC30vXj5aSQKmTJT907jQ7bIhD/ZXEfufV+9nJHJxomi8cH40
77Vu+4LSKC7uuksXjBpVkzvbx4Lu3+OC8fv0h+p4GinQs0gtt+Ui7Qsz6aZF
LclpdNiVJ8kF3VrUZ3NSLVl1ae54tQgXYpMqhzvYtUTceNu8cj4u1MO2XX9p
8o0kDN0Tc2DM4XGxizzF/EbuJ3VN+5XO4chqhmPRljqyX27TKpbDHPrtrx+X
sqwjCdaThy5az8H3vcjuy7fqSFjf2WU00zmkH/hbcs2POiLTonB07845iHnZ
5ftc/UH6dhi/4Vs0h6rGTiGdnHpyia0dlZ0wCxsF55HHHfWkrFNr6MWzWbAd
qRbGggaStLS1yvbeLNYLrkmNtGggisfmkzTnWdxRO2k2y24gjSf5ZiR2zELv
YcPLAu0momQTXiJeMYNCwyrNHaUtxNVyRzDnGwfM7dvK/gy3EL1AXUdTcLBG
IdokdUUrOUNrdQlJ5eAGx8VGzqWVeKVVtyUGc6D7dkWIkNxPoloZreWmy0H5
YoeOOp828tpIRPzaGzaqahfc/lerg5QHTq6XspkGp8BFQsm2g4xp+cpmGk9j
w5ufL+gPO8glkWh1VY1p3PNNSzvS30F0+cI/1y2aRuKzt4yIA53kIvFwzMia
gr2hmdQtTiepWEQ11s9NYig51uGgRRfJblu74v31CUy22gt8cuwiKXmX6hae
mwDfos0RO/27iMH1ieI9RhOQ+reAtv5dF7FP+axjLT2BPSr1m+fN7yYWny9k
78wax7O3gpM5ad1krDk+7mH3H+il2N+SX9BLJDxFLsduGkNU6qaI8Q99pN5w
QWB1GwOXF0sHsyv7iP2rpexTMQzs+pffg9veR55Zu5l/tGWAvr7eVFion3xe
srukZ3gIcvGeIquO9RNjs1HXtplBxL8o99w72U8k935QvyA2gJQ7lhZh6oMk
tqGKDIr2wvv3Xv2nRoPERrYsWpPWA5N9/2yLtB4k+d+W9prf7QFTgE8i8fYg
2XExV3lKsAeqNxK+5DUOErtrqwK7OF14f3VE/Zf7EPHlN3hl/rMTH8/eXqaa
wyBn7p9INx5uRotE+49eGoNE9o8f0TjSjEmohUV0MUi8MJ9syvsmqK75JbZA
bJhIPrXduc2rEXGtO4Xp54dJwLCjTxH1A3fNRjn3ljLJNF+fsftADV4LGOTt
Xs8kzYeluA9jaSjPiPaY0mGSwxuUhIotq8EVNx63cWASPq9fMicrKuBCSxzW
KmYS35D6Y4HPymCx71RnnyOL6M4YbDFfkQm38czoKH8WeRWxpWN92zuExS20
MnvOInwXDuVw/JKQZnpm5cI0Fum0GBpcrRgPGn9OS1E5z94lMq4YGIWBdNFn
V1t5fHpnZXp+GP7v/0j+///4P77AKTg=
      "]], LineBox[CompressedData["
1:eJxTTMoPSmViYGCwAGIQ3Tw9+EMR1wd7BihQuN/x8hQvgr/ObMuJpRIIvm3f
/eUNygh+lK1Zmpklgj959uNHC1MQfKXPfIdq8hH8DV5WC8OrEPwzPyfE805A
8Nki7O5U7Ebwq0RmXvUX+gjnZzOH3S+UQfBjPgm9nKyG4Ntf6Plz3QrBZ+lp
UE5IRvD7mDMKC7Yg+Is/mfFODPuEcM+FKy6XZn2G8wtPm6v1i32D85ecdxFN
O/EDzu9Zs30Fu89vOJ9x8hXPw6z/4HzxlauPL/FgcIDxrwFDj7WDEc6PXq3V
VvSICc63mc1ufVSVBc5/3L/zm0oLK5y/Rbn+teIdNjj/weJyIxZZDjj/koZr
w+cKTjhffparwom9XHD+tFWRPy3/cMP5l/c73dwRygvnZx+M45adwQfnJ82/
OW3PXAQ/qjZEKWYxgu9h6WU5ex2Cr7LRNF3yKIJ/ZwHPYdFPCL5Pw64qPh9+
ON85zpRtXSCCb2WzYaJvOIKv8WPpip4kBJ+lYOI1rkoEf298hhH7MgR/i+3j
vctWI/hrpOM83TYi+LOuBye07kHwy/zt+5gvI/i5urskF99A8FO4TZc63UPw
g45r7Wl8ieB7Ll3qrvgewXdoVrh84AuCb544Oy7hF4KvZy/26v9/BB8AV/7G
lw==
      "]]}, {}}, {{}, {}, {}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  Method->{
   "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None},
  PlotRange->{{-4, 4}, {-1.99999924218356, 0.}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.8227283366370473`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Integrate", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"an", "[", 
     RowBox[{
      RowBox[{"k", "-", "p"}], ",", "n"}], "]"}], "*", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"p", "^", "2"}], "-", "1"}], ")"}]}], ",", "p"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.8227281393416953`*^9, 3.8227281496727066`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{"n", " ", 
   RowBox[{"(", 
    RowBox[{"p", "-", 
     FractionBox[
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "4"}], "+", 
         RowBox[{"4", " ", 
          SuperscriptBox["k", "2"]}], "-", 
         SuperscriptBox["n", "2"]}], ")"}], " ", 
       RowBox[{"ArcTan", "[", 
        FractionBox[
         RowBox[{"2", " ", 
          RowBox[{"(", 
           RowBox[{"k", "-", "p"}], ")"}]}], "n"], "]"}]}], 
      RowBox[{"2", " ", "n"}]], "+", 
     RowBox[{"k", " ", 
      RowBox[{"Log", "[", 
       RowBox[{
        SuperscriptBox["n", "2"], "+", 
        RowBox[{"4", " ", 
         SuperscriptBox[
          RowBox[{"(", 
           RowBox[{"k", "-", "p"}], ")"}], "2"]}]}], "]"}]}]}], ")"}]}], 
  RowBox[{"2", " ", "\[Pi]"}]]], "Output",
 CellChangeTimes->{3.822728150677638*^9}]
}, Open  ]]
},
WindowSize->{1536, 781},
WindowMargins->{{-8, Automatic}, {Automatic, -8}},
Magnification:>2. Inherited,
FrontEndVersion->"10.1 for Microsoft Windows (64-bit) (March 23, 2015)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 399, 12, 66, "Input"],
Cell[982, 36, 496, 18, 66, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[1515, 59, 819, 28, 69, "Input"],
Cell[2337, 89, 735, 27, 100, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3109, 121, 135, 2, 60, "Input"],
Cell[3247, 125, 580, 21, 99, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3864, 151, 600, 19, 60, "Input"],
Cell[4467, 172, 514, 18, 115, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5018, 195, 97, 1, 60, "Input"],
Cell[5118, 198, 494, 18, 115, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5649, 221, 176, 4, 60, "Input"],
Cell[5828, 227, 593, 20, 142, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6458, 252, 215, 5, 60, "Input"],
Cell[6676, 259, 591, 21, 99, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7304, 285, 153, 3, 60, "Input"],
Cell[7460, 290, 437, 15, 91, "Output"]
}, Open  ]],
Cell[7912, 308, 526, 17, 60, "Input"],
Cell[8441, 327, 253, 7, 60, "Input"],
Cell[CellGroupData[{
Cell[8719, 338, 558, 16, 60, "Input"],
Cell[9280, 356, 1301, 43, 184, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10618, 404, 184, 3, 60, "Input"],
Cell[10805, 409, 1287, 43, 143, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12129, 457, 268, 7, 60, "Input"],
Cell[12400, 466, 916, 30, 141, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13353, 501, 121, 2, 60, "Input"],
Cell[13477, 505, 98, 1, 59, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13612, 511, 173, 4, 60, "Input"],
Cell[13788, 517, 508, 17, 113, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14333, 539, 332, 10, 60, "Input"],
Cell[14668, 551, 1220, 42, 168, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15925, 598, 129, 2, 60, "Input"],
Cell[16057, 602, 1220, 42, 168, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17314, 649, 135, 2, 60, "Input"],
Cell[17452, 653, 70, 1, 59, "Output"]
}, Open  ]],
Cell[17537, 657, 353, 11, 60, "Input"],
Cell[17893, 670, 285, 9, 60, "Input"],
Cell[CellGroupData[{
Cell[18203, 683, 390, 11, 60, "Input"],
Cell[18596, 696, 3751, 74, 476, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22384, 775, 201, 5, 60, "Input"],
Cell[22588, 782, 583, 19, 134, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23208, 806, 250, 7, 60, "Input"],
Cell[23461, 815, 845, 26, 90, "Output"]
}, Open  ]],
Cell[24321, 844, 678, 22, 141, "Input"],
Cell[CellGroupData[{
Cell[25024, 870, 148, 3, 60, "Input"],
Cell[25175, 875, 387, 12, 109, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25599, 892, 346, 10, 60, "Input"],
Cell[25948, 904, 398, 11, 92, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26383, 920, 124, 2, 60, "Input"],
Cell[26510, 924, 131, 2, 59, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26678, 931, 760, 23, 99, "Input"],
Cell[27441, 956, 212, 3, 59, "Output"]
}, Open  ]],
Cell[27668, 962, 454, 14, 60, "Input"],
Cell[CellGroupData[{
Cell[28147, 980, 744, 21, 60, "Input"],
Cell[28894, 1003, 13872, 241, 540, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[42803, 1249, 405, 12, 60, "Input"],
Cell[43211, 1263, 997, 30, 99, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[44245, 1298, 275, 8, 60, "Input"],
Cell[44523, 1308, 394, 12, 90, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[44954, 1325, 488, 15, 60, "Input"],
Cell[45445, 1342, 77, 1, 59, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[45559, 1348, 524, 15, 60, "Input"],
Cell[46086, 1365, 868, 27, 99, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[46991, 1397, 135, 2, 60, "Input"],
Cell[47129, 1401, 1041, 34, 111, "Output"]
}, Open  ]],
Cell[48185, 1438, 484, 15, 60, "Input"],
Cell[CellGroupData[{
Cell[48694, 1457, 145, 3, 60, "Input"],
Cell[48842, 1462, 302, 9, 136, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[49181, 1476, 491, 13, 60, "Input"],
Cell[49675, 1491, 1516, 45, 192, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[51228, 1541, 131, 2, 60, "Input"],
Cell[51362, 1545, 1443, 45, 192, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[52842, 1595, 1031, 33, 99, "Input"],
Cell[53876, 1630, 914, 30, 98, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[54827, 1665, 321, 10, 60, "Input"],
Cell[55151, 1677, 1142, 38, 121, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[56330, 1720, 181, 4, 60, "Input"],
Cell[56514, 1726, 378, 12, 89, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[56929, 1743, 353, 10, 60, "Input"],
Cell[57285, 1755, 124, 2, 59, "Output"]
}, Open  ]],
Cell[57424, 1760, 1934, 60, 175, "Input"],
Cell[CellGroupData[{
Cell[59383, 1824, 254, 6, 60, "Input"],
Cell[59640, 1832, 194, 3, 59, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[59871, 1840, 71, 1, 60, "Input"],
Cell[59945, 1843, 1122, 38, 121, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[61104, 1886, 170, 4, 60, "Input"],
Cell[61277, 1892, 1259, 43, 124, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[62573, 1940, 181, 4, 60, "Input"],
Cell[62757, 1946, 1868, 60, 211, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[64662, 2011, 122, 2, 60, "Input"],
Cell[64787, 2015, 159, 3, 59, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[64983, 2023, 71, 1, 60, "Input"],
Cell[65057, 2026, 1120, 38, 121, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[66214, 2069, 170, 4, 60, "Input"],
Cell[66387, 2075, 1243, 41, 121, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[67667, 2121, 391, 11, 60, "Input"],
Cell[68061, 2134, 123, 2, 59, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[68221, 2141, 123, 2, 60, "Input"],
Cell[68347, 2145, 1355, 43, 121, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[69739, 2193, 216, 6, 60, "Input"],
Cell[69958, 2201, 393, 13, 99, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[70388, 2219, 144, 3, 60, "Input"],
Cell[70535, 2224, 255, 8, 111, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[70827, 2237, 69, 1, 60, "Input"],
Cell[70899, 2240, 1122, 38, 121, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[72058, 2283, 154, 3, 60, "Input"],
Cell[72215, 2288, 471, 10, 46, "Message"],
Cell[72689, 2300, 82, 1, 59, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[72808, 2306, 71, 1, 60, "Input"],
Cell[72882, 2309, 1122, 38, 121, "Output"]
}, Open  ]],
Cell[74019, 2350, 71, 1, 60, "Input"],
Cell[74093, 2353, 1399, 44, 191, "Input"],
Cell[75495, 2399, 271, 9, 60, "Input"],
Cell[CellGroupData[{
Cell[75791, 2412, 218, 4, 60, "Input"],
Cell[76012, 2418, 187, 4, 59, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[76236, 2427, 170, 4, 60, "Input"],
Cell[76409, 2433, 117, 3, 91, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[76563, 2441, 122, 2, 60, "Input"],
Cell[76688, 2445, 91, 1, 59, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[76816, 2451, 220, 6, 60, "Input"],
Cell[77039, 2459, 351, 12, 111, "Output"]
}, Open  ]],
Cell[77405, 2474, 2255, 70, 246, "Input"],
Cell[CellGroupData[{
Cell[79685, 2548, 179, 4, 60, "Input"],
Cell[79867, 2554, 2094, 68, 211, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[81998, 2627, 234, 5, 60, "Input"],
Cell[82235, 2634, 141, 2, 59, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[82413, 2641, 355, 10, 60, "Input"],
Cell[82771, 2653, 144, 2, 59, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[82952, 2660, 897, 28, 213, "Input"],
Cell[CellGroupData[{
Cell[83874, 2692, 189, 2, 43, "Print"],
Cell[84066, 2696, 188, 2, 43, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[84303, 2704, 373, 9, 60, "Input"],
Cell[84679, 2715, 2147, 70, 236, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[86863, 2790, 279, 7, 60, "Input"],
Cell[87145, 2799, 4446, 85, 508, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[91628, 2889, 519, 14, 60, "Input"],
Cell[92150, 2905, 2430, 70, 196, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[94617, 2980, 273, 7, 60, "Input"],
Cell[94893, 2989, 691, 22, 152, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[95621, 3016, 664, 19, 60, "Input"],
Cell[96288, 3037, 148, 2, 84, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[96473, 3044, 418, 13, 60, "Input"],
Cell[96894, 3059, 327, 11, 126, "Output"]
}, Open  ]],
Cell[97236, 3073, 120, 2, 60, "Input"],
Cell[CellGroupData[{
Cell[97381, 3079, 473, 15, 60, "Input"],
Cell[97857, 3096, 7088, 130, 501, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[104982, 3231, 344, 11, 60, "Input"],
Cell[105329, 3244, 861, 28, 166, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
